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Abstract 

We present results about large deviations and laws of large numbers for various polymer 
related quantities. 

In a completely general setting and strictly positive temperature, we present results 
about large deviations for directed polymers in random environment. We prove quenched 
large deviations (and compute the rate functions explicitly) for the exit point of the 
polymer chain and the polymer chain itself. 

We also prove existence of the upper tail large deviation rate function for the loga- 
rithm of the partition function. In the case where the environment weights have certain 
log-gamma distributions the computations are tractable and allow us to compute the 
rate function explicitly. 

At zero temperature, the polymer model is now called a last passage model. With 
a particular choice of random weights, the last passage model has an equivalent rep- 
resentation as a particle system called Totally Asymmetric Simple Exclusion Process 
(TASEP). We prove a hydrodynamic limit for the macroscopic particle density and cur- 
rent for TASEP with spatially inhomogeneous jump rates given by a speed function that 
may admit discontinuities. The limiting density profiles are described with a variational 
formula. This formula enables us to compute explicit density profiles even though we 
have no information about the invariant distributions of the process. In the case of a 
two-phase fiux for which a suitable p.d.e. theory has been developed we also observe 
that the limit profiles are entropy solutions of the corresponding scalar conservation law 
with a discontinuous speed function. 
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Chapter 1 
Introduction 

1.1 Polymers at finite temperature and corner growth 
models 

We begin by presenting the two main models that are discussed in this dissertation. 
After the two models are introduced we offer a connection between the two of them 
(namely one can be viewed as a limiting case of the other) . The two remaining sections 
of the chapter can be viewed as an informal introduction to the material that follows: 
Some basic definitions, discussion on classical results and an idea of the kind of questions 
we are asking. At the end of the chapter we describe the organization of the thesis. 

1.1.1 General polymer models 

A directed polymer in random environment is a random walk path that interacts with 
a random environment. The polymer chains live in x Z_|_, where the last coordinate 
denotes time. The space of environments is denoted by = {a;(u, n) : u e Z*^, n e Z+} 
and is equipped with a probability measure P, so that under P the random variables 
a;(u, n) are i.i.d. for all u, n. 

The two models under consideration are directed polymers with free endpoints and 
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directed polymers with constrained endpoints. Here, directed means that the last coordi- 
nate is always increasing by 1 at each time step. This allows for the polymer chain in 
d + 1 dimensions to be viewed as the path of a d-dimensional nearest neighbor random 
walk. 

We assume that at time t = 0, the starting point of the polymer chain (the random 
walk) is anchored at e Z'^. For each m e N, define the set of possible endpoints for the 
polymer to be £{m). For a fixed (u, m) in S{m), the set of all polymer chains starting 
from (0, 0) — Xq and ending at (u, m) — is 

7^(U, m) = {xo,m ■Xo,m = (^^O, 2^1, ■■■ , Xm), Xk - Xk-1 = (±6^, 1) (1.1) 

for some 1 < i < d, for all 1 < /c < m, (0, 0) = Xq, (u, m) — x^}, 

where Ci,! < i < d is the standard basis of M."^. 

The point-to-point partition function is defined by 

m 

and the total (or the point-to-line) partition function can be defined by 

m 

E E n^"^^^^^ (1-3) 

The parameter P is what is known in the literature as the inverse temperature and 
is assumed without loss of generality to be positive. 

Under a fixed environment ou, the polymer chain xo,m is selected according to the 
quenched probability measures 

m 

Qu:i{^o,m} = {ZP{n,m))-'l[e^-(^^\ (1.4) 

i=i 
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and 

m 

Q^J{^o,m] = {Z^{m)r'Y[e^^^^^\ (1.5) 

respectively for each of the models described above. 

Let us momentarily restrict our attention to dimension 1 + 1. The polymer chain 
starts by being anchored at (0, 0) and under a fixed environment, at time n the chain is 
chosen according to the measures Q'^'^ . The chain lives inside the cone {u = (mi,M2) G 
Z X Z_|_ : U2 > \ui\}. We rotate the picture clockwise by 45 degrees. The polymer chain 
now becomes a path xo^rn = {(0, 0) = Xq, Xi, • • • , Xm} in the first quadrant where the 
differences Xk — Xk-i are a standard basis vector. Such a path we call an up-right path. 

The advantage of this viewpoint resides in the point-to-point model: For any vector 
u e we specify as an endpoint, it is guaranteed that exponentially many paths 
start at (0, 0) and end at u, i.e. u is an admissible endpoint. Unfortunately we lose 
information about the time (now the time axis is the main diagonal) but this does not 
affect the analysis and the limits of the main theorems. 

Prom this point onwards, all models in this dissertation are about paths that live 
in the first quadrant and are up-right paths. We denote the set of up- right paths from 
(0, 0) to (u) by n(u). More precise details are available in the following chapters. 

1.1.2 Corner growth model 

The corner growth model in two dimensions is a model of a randomly growing cluster 
that over time covers larger and larger portions of the first quadrant Z^. At the outset, 
each coordinate is given a weight uj{i,j). In the language of the previous section, 
we have now a fixed environment and we want to study the evolution of the random 
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Figure 1: The shaded squares represent sites that are occupied at this particular time. 
The thick hne represents the boundary of the cluster and the squares marked with X 
are the " growth sites " . The time it takes each growth site to be infected (shaded) is 
1^0,451^1,2,1^2,1,1^3,0 respectively. These times start elapsing as soon as both the left and 
lower neighbors become shaded. 

cluster according to the following rules. 

The general rule is that uj{i, j) is the time it takes for the random evolution to occupy 
site This can only be done only after its two neighbors to the left and below are 

either occupied or lie outside Z^. At the boundaries the rule is that point (0, 0) needs no 
occupied neighbors to start, points (1, j) on the left boundary wait only for the neighbor 
below to be occupied, and points (i, 1) on the bottom boundary wait only for the left 
neighbor to be occupied (see Fig. 1). 

The quantity of interest is that of the last passage time. Given a site the last 
passage time, T{i,j), is the time when the site becomes part of the evolving cluster. 
Using the evolution described above, it can be recursively computed (with appropriate 
boundary conditions) to be 

T(^,j) = T{i - l,j) - 1) +0;^,,. (1.6) 

Equation (1.6) says that becomes part of the cluster only after both [i — 1, j) and 
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— 1) are part of the cluster, and after that happens, time uuij elapses. Assume 
without loss that T(0,0)=a;(0,0)=0. An easy induction argument then yields that 

i+j 



In the case where the ou weights are continuous and an endpoint u = (i, j) is specified, 
there is a unique path a;o,i+j for P- a.e. u that attains the last passage time. Wc call that 
the maximal path and denote it by x^f^j- In this respect, we can define a (degenerate) 
quenched measure on the paths. 



If the u weights have discrete distributions, is possible that there are more than one 
maximal path. When that happens, the quenched probability measure on the paths is 
the uniform measure on maximal paths. 

1.1.3 Connection of the two models 

Now we are ready to justify the title of this dissertation. The connection between the 
polymer models and the last passage time models comes via the parameter (3. As we 
let (3 tend to oo, under a fixed environment, the quenched probability measure Q'^'^ 
converges weakly to a delta mass, given by (1.8). This becomes precise in the context 
of the next proposition. 

Proposition 1.1. Let u e Z^, u = and let m — i + j. Fix an environment ou. 

Then, the probability measures defined by (1.4) converge weakly 




(1.7) 




(1.8) 




as P ^ oo. 



(1.9) 
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Proof. Let / a function on the space n(u) of up-right paths. Then 



m 



a:o,men(u) j=i 

m 



a;o,men(u) j=i 

m 

fK^)Z^n,mr'lle^-K-) (1.10) 



X 



(l+ E fMe^'^^'^-^^^^-^^^'^A (1.11) 

\ a;o,mSn(u) / 

As /3 — >■ oo, (1.11) tends to 1 because all the exponents are negative. We are going to 
show that Zl^(u, m)~^ YYjLi e'^'^^^?"'''-' tend to 1 as ^ — >■ oo. Then the proposition follows 
from (1.10) and (1.11). 

Prom the definitions, Z^{u, m)"^ YYjLi e^^^^i"^""^ < 1. For a lower bound, observe that 
for any S > and (3 sufficiently large 

en(u) 

<l + 5. □ 
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1.2 Large deviations 

The theory of large deviations is concerned with the study of rare (improbable) events. 
Assume you have a sequence {X„}jjgN of real random variables on Q. A natural question 
that arises is to compute the limiting probability limn-^^¥{Xn G A} where A is a 
fixed measurable set on R. When the events An — {Xn G A} are "rare" the limiting 
probability can be 0. This information, while not particularly helpful, leads to a deeper 
question: How fast does it go to 0? For example if the probabilities are now summable, 
one can apply the Borcl-Cantelli lemma. 

In a more general setting, instead of a sequence of random variables, one can use a 
sequence of probability measures //„ on a Polish measurable space X (in the example 
above, /x„ = P{X„ e •} and A* = M). The measures /x„ live in the space of probability 
measures on X, A4i{X) and we say that they satisfy a large deviation principle if the 
following definition holds. 

Definition 1.2. Let / : A" i->- [0, oo] be a lower semicontinuous function and r„ — >■ oo 
a sequence of positive constants. A sequence of probability measures {/x„}„ C A^i(A') is 
said to satisfy a large deviation principle with rate function I and normalization rn if 

the following inequalities hold for all closed F C X and all open G <^ X: 



When the sets {/ < c} are compact for all c e R, we say / is a tight rate function. 
It is of interest to find explicit rate functions, since they offer an exact measurement of 




lim — log//„(F) < — infJ; 



(1.12) 



lim — log/in(G) > — inf /; 



(1.13) 
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the rare event. For example, insurance companies can use that information to decide on 
a fair premium for the customer. 

1.3 Interacting particle systems and hydrodynamic 
limits 

In full generality, interacting particle systems consist of finitely or infinitely many par- 
ticles that evolve in space and time according to given transition probabilities or rates, 
with some interaction rules imposed on the particles. A particle system that has been 
extensively studied, is the Totally Asymmetric Simple Exclusion Process (TASEP). It is 
intimately connected with the last passage time, assuming exponential weights. 

Assume that particles occupy integer sites. There is at most one particle at each site 
(Simple). Each particle attempts to move one unit to the right (Totally Asymmetric) 
with rate 1 . The jump is suppressed with probability 1 if the target site is occupied 
(Exclusion). 

One-dimensional TASEP can be constructed graphically by assigning independent 
mean 1 Poisson processes (called clocks) on each integer site. The vertical direction now 
becomes time. As time progresses, a particle on site i attempts to jump at the Poisson 
event times of the site it occupies, and the jump happens with probability 1 as long 
as the exclusion rule is satisfied. With probabihty 1, two adjacent Poisson processes 
cannot have simultaneous events, and for every time t there are infinitely many Poisson 
processes with no events before t, so one can study the temporal evolution of the system 
up to time i in a rectangle [— M, N] x [0, t] around the origin. 



9 



The coupling with the corner growth model follows if we run a TASEP starting from 
step initial conditions: particles start by occupying only the negative integers and are 
labeled so that particle i starts on site —i. Then the last passage time r(m, n) is equal 
in distribution to the time it takes the n-th particle to reach site m — n. 

Consider a sequence of exclusion processes rp = (?7"(t) : i G Z, t G M+) indexed by 
n E N. For each i and fixed t, r)f{t) = 1 if and only if there is a particle present at 
site i, at time t. These processes are constructed on a common probability space that 
supports the initial configurations {77" (0)} and the Poisson clocks of each process. The 
clocks of process r^" are assumed to be independent of its initial state 77" (0). 

Starting from arbitrary particle initial conditions and for fixed time t define the 
sequence of occupation measures 

1 ^"'^ 1 
li^{[a,b]) = — ^ V?{'>T't) = particles in [na,nb] at time nt}. (1-14) 

i=\na\+l 

Under some regularity assumptions on the initial conditions, it is known (e.g. [27]) 
that this sequence of measures converges weakly for all t 



where p(x, t) (called the particle density function) is the unique entropy solution to the 
scalar conservation law 



The initial conditions of the particles correspond to the initial conditions required for 
uniqueness of the weak solution in (4.19) and F is the particle flux of TASEP, given by 
F(x, p) — p{l — p). Results of this type go by the name of hydrodynamic limit and there 
are many known generalizations. 




(1.15) 



dt{p{x,t)) + d,F{x,p{x,t)) = 0. 



(1.16) 
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1.4 Motivation 

Upper tail large deviations for the last passage time have been computed in the case of 
geometric and exponential weights (see [17, 26]). In [26], the rate function was computed 
via the height function and information about equilibrium distributions for TASEP 
particles. The equilibrium distributions for TASEP is a result of Burke's theorem for 
M/M/1 queues and they can be interpreted as appropriate boundary weights in the 
corner growth model. This (Burke) property is "transferable" to the log-gamma polymer 
model. The model was introduced in [29]. The model is a 1 + 1 dimensional polymer 
model at temperature (3 — 1 where for a fixed > 0, the weights 

^{hj) — ~log^j) ~ Gamma(//) 

for (i,j) e N2. 

1.5 Organization 

We start (Chapter 2) by describing general polymer models. In Chapter 2 we show 
some general facts about the partition function and proceed by showing quenched large 
deviation results for the polymer chains under the quenched measures (1.4), (1.5). In 
a completely general setting we show some quantitative properties of the rate functions 
(existence, continuity, large (3 behavior). 

In Chapter 3 we restrict to a specific 1 + 1 dimensional model, the log-gamma model 
that satisfies a certain property that allows for explicit computations. We present large 
deviation results about the logarithm of the partition functions. 

The following two chapters (Chapters 4 and 5) are concerned about hydrodynamic 



11 



limits of exclusion processes and last passage time where the weights u are now ex- 
ponential but with different parameters (so they are not identically distributed). The 
parameters are decided using a (possibly discontinuous) function c{x). The connection 
with the particle process TASEP leads to a further connection with scalar conservation 
laws with discontinuous coefficients. 
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Chapter 2 

Generalities about Polymer Models 



2.1 Introduction 

2.1.1 Directed polymers with constrained endpoint in a rect- 
angle 

For u G define the set of directed polymer chains from v to u with ||u — v||i = m 

nv(u) = [x. = {v = Xq, xi, . . . , Xm = n} : Xk e Z'^, 

Xjfc+i — Xk — ei for some 1 < i < d}, (2.1) 

where is the i-th standard basis vector. 

The point-to-point partition function is in this case defined by 

m 

^v,u= E n^''"^'''^= E e'^^r=i-(-.). (2.2) 
a;.env(u)i=i a;.env(u) 

In the special case where v = we omit the index v from the above notation: The 

partition function is denoted by and the set of polymer chains by n(u). Observe 

that in the definitions given so far, the weight at the starting point is ignored. When 

that is not the case, we denote 
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Under a fixed environment a;, fixed endpoint u and fixed inverse temperature ^, the 
quenched probability measure Qu'^ on paths with constrained endpoints, is defined by 

m 

Q-J{x.) = {Z^)-'l[e^-^-^\ (2.3) 

2.1.2 Directed polymers with constrained endpoint in a rect- 
angle 

In this variation we fix the number of time steps m. Define the set of aU admissible 
polymer chains starting from to be 

IltotM = {xo,m ■ Xo^m = {xq, Xi, , ... Xm} , Xk+1 - Xfc = 6^ (2.4) 

for some 1 < i < d, for all 1 < /c < m} 

where Si,! < i < d is the standard basis of M*^. 

For each m > 0, the total partition function is defined by 

ZT'- E ^u- (2.5) 

u6Z^:||u||i=m 

The corresponding quenched probability measure Q'^^ on paths in ntot(^), is 

m 

Q^JM = {z^"'T'll^^''''^'^- (2-6) 

i=i 

Remcirk 2.1. // the rectangle has dimension 2, then both models describe the classical 
directed random polymers in random environment in dimension 1+1 described in the 
introduction (point-to-point and free endpoint respectively), where the picture is rotated 
by 45 degrees, so the polymer lives in the first quadrant. 
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2.1.3 Known results 

Concentration inequalities. A problem in the area of random polymers in random 
environment is about the fluctuations (in particular the fluctuation exponent) of log 
that is conjectured to be 2/3 in the physics literature. Rigorous results about upper 
and lower bounds for fluctuations exponents for specific polymer models can be found 
in [8, 12, 20, 22, 33, 34]. The only two cases where the conjectured value of 2/3 is in 
fact verified is the log-gamma polymer [29] and the brownian polymer [30] . 

In absence of exact variance bounds, information can be derived from concentration 
inequalities. In [9] an exponential of order n^/^ concentration result has been obtained 
for the partition function in a Gaussian environment and later this had been generalized 
for all weights under certain exponential moment assumptions in [11]. 

The sharpest concentration inequalities are exponential of order n, proven in [19, 32] 
and a certain version used here about the point-to-point free energy in [13]. 

Large deviations. The exponential order n concentration inequality for the parti- 
tion function is the correct one for the upper tail large deviations and is the correct 
normalization in order to get a non-trivial upper tail rate function. 

For the lower tail the behavior depends on the distribution of the weights. In [7] 
three different normalization regimes for the lower tails are shown, depending on the 
distribution of the weights. It is also mentioned that a normalization is true if the 
weights are bounded. This is proved in the case where the weights are Gaussian or 
bounded in [9] where upper and lower normalizations are proven. 

In the last part of the chapter we prove quenched large deviations for the exit point 
of the polymer (in the free endpoint case) and quenched large deviations for the polymer 
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path in both the free endpoint and constrained endpoint. The results are described in 
further detail in Section 4.2. 

Notation and conventions. Throughout we use N for positive integers, while Z+ 
denotes the set of non-negative integer. Similarly, R+ denotes the set of non-negative 
real numbers and M.'^ is the set of all vectors with non-negative coordinates. All vectors 
{vi,V2t . . iVd) e M'' are denoted by boldface notation v = {vi,V2t ■ ■ iVd). The ordering 
V < u means vi < ui,V2 < U2,..-,Vd < ua- The d-dimensional vector with entries 
equal to 1 is denoted by 1 = (1, 1, . . . , 1) and correspondingly, = (0, 0, . . . , 0). For the 
purposes of this dissertation we define for y = (yi, y2, ■ ■ ■ , Vd) £ K*^, n eN the floor of a 
vector [nyj = ([nyij, L^?/2j, • • • , [nyd])- 



Some technicalities before stating the general results. We need a technical assumption 
on the environment to. Henceforth we assume 

Assumption 2.2. There exists some ^ > that depends on the distributions of the 

weights uj such that 



2.2 Definitions and Statement of Results 



(2.7) 



This guarantees the existence of a large deviation rate function defined by 



n 



I[r) = — lim lim n logP 



r 



— s,r 




(2.8) 



i=l 
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All results that follow are valid for ^ < ^. In order to make the proofs cleaner we also 
assume that for all u e Z^, 



We start by proving some qualitative properties for the limits of the rectangle par- 
tition functions and we summarize them in the following two propositions. The super- 
script 13 is considered fixed and henceforth omitted until the end of the chapter where 
it becomes relevant to take limits as /3 ^ oo. 

Proposition 2.3. Let y = 1/2, ... , yd) G W^, n e N and let Zy^y^ he defined by (2.2). 
Then the limit 



Furthermore, there exists an event Qq C Q of full probability such that the a.s. conver- 
gence happens simultaneously for all y e R^. The limiting value viewed as a function 
of y satisfies the following properties: 

(a) p(-) is concave and continuous on M!|_. 

(b) p(cy) = cp(y) for c> 0. 

(c) For all y such that \\y\\i = 1, p(y) < p{d~^l). 

Proposition 2.4. Let n G N, t > and let Zf"*, be defined by (2.5). Then the limit 



P{c^(u) > 0} > 0. 



(2.9) 



hm n ^ log Zyny\ = p{y) P - a.s. 



(2.10) 



lim n-^ log = p'°\t) P - a.s. 



(2.11) 



Furthermore, 



sup p(y) = p{d ^tl) 



(2.12) 



y:||y||i=< 
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These propositions suffice to guarantee the following general existence theorems. 
Then, some quantitative properties of the rate functions follow. 

Theorem 2.5. Under assumption 2.2, the following large deviations principles hold. 

1. For t > 0, u e \ {0} and r e R there exists a nonnegative function J that 
satisfies 

Ju(r) = - lim n-^ logP{log > nr}. (2.13) 

n— >-oo 

J is convex in the variable (u,r) and equals on the set r < p(u). The rate 
function is continuous in (u, r) on the set int{(u, r) : u e \ {0}, r e R, Ju{f) < 
oo} C R^ X R. 

2. For t > and r e R, there exists a nonnegative function J that satisfies 

Jt{r) = - lim n-MogP{logZ[°*| > nr}. (2.14) 

J is convex in the variable {t,r) and equals on the set r < p^°^{t). The rate 
function is continuous in {t,r) on the set mt{{t,r) : Jt{r) < oo}. 

Next, we show some quantitative properties of the rate functions described in the 
following three propositions. 

The first one is about the behavior of the rate function Ju{r) at the lower di- 
mensional boundaries of R^. For a vector u = {ui,...,Ud) e R^ we denote by 
Ui,A; = {ui, . . . , Uk, 0, . . . , 0) the projection onto the k-dimensional facet of the first quad- 
rant. 

It is possible that u = so we need a definition for the rate function at 0. If we just 
change u with in (2.13) the rate function becomes trivial (takes values and oo only). 
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However, 0, is the macroscopic endpoint of the partition function. The definition of Jq 
should cover for example the rate function of a single random variable. It is consistent 
(in the sense that limu^o Jui^^) = •^(^) described in the following proposition) to define 
it the following way: 

{0, r < 0, 
(2.15) 
Xoor, r > 

where e (0, oo] is the maximal slope of the one-sided Cramer rate function on the 
right of the zero for sums of i.i.d. ou weights and is given by Xoo — linir->oo r~^/(r). 

Proposition 2.6 (Continuity at the boundaries). Let u e and fix r e R. Assume 
that Jiiuiiieil'") < oo in a neighborhood of ||u||iei. (This is the one sided Cramer rate 
function.) Then, for any k > 1 and sequence {u^'^^j^gN where u*^"*^ < u, lim^_^oo u^'"-' = 

lim J^(m){r) = Ju^i^ir). (2.16) 

m— >-oo 

In the special case where Jo{r) < oo, the function J{.){r) is continuous everywhere on 
R^. 

Proposition 2.7 (Unique zero). 

1. The rate function Ji{r) given by (2.14) is strictly positive for r > p*°*(l). 

2. The rate function Ju('^) given by (2.13) is strictly positive for r > p(u). 



As stated in the introduction, the polymer model starts behaving hke the last passage 
time model when (3 ^ oo. For a fixed u, define the last passage time as in (1.7). 
Assumption 2.2 along with the subadditive ergodic theorem, give the existence of a 
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finite limiting last passage time constant 



lim n~^T{[nu\) ^ip{u), 



with arguments identical to those used in the proof of Proposition 2.3. In turn, one can 
show the existence of an upper tail large deviation rate function for the last passage 
time, following the steps of the proof of Theorem 2.5. The rate function is given by 



The rate functions also start behaving similarly as described in the following propo- 
sition. 

Proposition 2.8 (Large /3 behavior). Let u e and assume I^{r) is left continuous 
at r eM.. Then 



We also prove quenched large deviations for the exit point for the free endpoint 
models, and the path of a polymer chain for both models. 

Theorem 2.9 (Exit point LDP). Let u = {ui,U2, ■ ■ ■ ,Ud) G with ||u||i = 1 and let 
y — (iti, . . . , iid_i) e R^"-^. For n e'N let [nu] — {[ny\,n — \\ [ny\ ||i) and denote by Xn 
the last point of the polymer chain xo,n- Then 




lim n ^ log P{T( [nuj ) > nr}. 



lim J^i/Sr) = /~(r). 



(2.17) 



p{d-'l)-p{u) 



lim n ^ log(55^{x 



n 



[nu]}. 



(2.18) 



This readily leads to the following path large deviations. 
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Theorem 2.10 (Path LDP). Let 7 : [0,1] — > be a coordinatewise nondecreasing 
Lipschitz curve with 7(0) = 0. For e > let Me{l) denote a uniform e -neighborhood of 
7. The following large deviation principles hold: 

1. (Constrained endpoint.) Let u in and let 7(1) = u. Then, 

p{n) - I p{j'{t))dt = - Urn Urn n'HogQr^^Ax. e nK{l)} (2.19) 

Jq £->0n— >-oo l j <- 

2. (Free endpoint.) Let ||7(l)||i = 1 • Then, 

p{d-H)- [ p{^'{t))dt^-limlimn-HogQ'^{xo,nenKil)} (2.20) 

Remcirk 2.11. The theorem is true for any curve 7 that admits a Lipschitz parametriza- 
tion. It is easy to check that the rate functions in (2.10) are independent of any 
parametrization (f) : [0, 1] — >■ [0, 1]. It follows from the 1- homogeneity of p and a change 
of variables in the integral. 



2.3 Preliminaries 

In this section we record resuhs that are used throughout. In particular, basic properties 
of the partition functions; their laws of large numbers and the proof of Propositions 2.3, 
2.4. For the continuity of the partition functions at the boundary and for the unique 
zeros of the rate functions, we need a concentration inequality which we state first. At 
the end of the section we prove an auxiliary lemma about upper-tail large deviations of 
a sum of independent random variables and conclude the section with a basic fact that 
we use at various instances without alerting the reader. 
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Lemma 2.12 (Proposition 3.2.1 [13]). Under assumption 2.2, there exist constants 
Ci, C2 e (0, oo) such that for all e > sufficiently small 

P{|logZ„ -ElogZ^I > ne} < 2exp{-ci£^n} (2.21) 

and 

P{|logZL„„j -Elog^L^ujI > ne} < 2exp{-C2{n)e^n} (2.22) 



Proof. The proof is identical with the one in [13], p. 26, tailored to these particular 
rectangle models. □ 

For y e Z^, define the coordinate shift operator Ty so that 

Zu,V °Ty — Zu+y,v+y (2.23) 

and observe that the random variables log Zu,v are super additive, satisfying 

log > log + log (2.24) 

= log + log Zv O Tu. 

Proof of Proposition 2.3. Assume without loss of generality that K{uj{0)) > 0. This 
builds the monotonicity of the limiting free energy that will simplify the proof and makes 
the limiting free energy a positive function. The proof for concavity and homogeneity 
when y e N'^ and c e N follows from the subbadditive ergo die theorem. For the 
monotonicity of the limit, let u < v e and for given n fix a path 7r„ from nu to nv. 
Superadditivity gives 

log > log Z^^ + /3J2 ^(^i)- (2-25) 



22 

Taking a limit along a suitable subsequence after dividing by n gives that p(v) > p(u) + 
||v — u||iE(c<;(0)), hence p is monotone on integer vectors. 

For rational y G Q*^, find a positive integer k so that ky e N*^ and define 

p(y) = k-^ky)- (2.26) 

Observe that homogeneity for integer k gives that the value in (2.26) is independent of 
the choice of k. Homogeneity and concavity extend now to rational c > and rational 
y as follows. 

Let n e N and write it as n = Mk + r with r e (0, 1, . . . , A; — 1}. Then 

Mky < [Mky + ryj = [nyj < (M + l)ky. (2.27) 
Hence, by the superadditivity 

-^(M+l)fey ^'^[nyJ,(M+l)fey^ > -^[nyj > -^Mfey -^Mifcy,[nyJ ■ (2.28) 

We show the upper bound. The lower one follows in the same manner. Observe that 
||(M + l)A;y — ny||i < A;||y||i, so as n — >■ oo, n~^log Z^nyi,{M+i)ky — >■ almost surely. 
This, along with (2.28), gives 

^ I ^ 

p{y) = lim — — - log Z(„_^+fe)y > lim -logZ^nyj. (2.29) 

Monotonicity follows as in the case with integer coordinates. 
To extend to all vectors y in M^, define 

p{y) = sup {p(x„) : y > x^ e Q''} (2.30) 

With this definition, homogeneity follows immediately for rational c. For any c > pick 
ci,C2 e Q, so that ci < c < 02- Then 

cip(y) = p(ciy) < p{cy) < p{c2y) = C2p(y). (2.31) 
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Then consider sequences c and \ c to get homogeneity for all c e R+. 

Superadditivity follows directly from (2.30) and concavity follows from homogeneity 
and superadditivity. This in turn implies continuity of p(-) in the interior of R^. 

Finally, for the most general form of the limit for y e R^, pick rational vectors 
Ui < y < U2 so that for a given s > 0, \\ui — y||i < £■ By superadditivity, 

log Zlnui] + log ^[nuiJ,LnyJ < log ^[nyj < log ^LnuaJ " log ^[„yj,[nu2j • (2-32) 

We first bound log Ziny],[nu2i from below using a path 7r„ with segments on the 
boundary of the rectangle, parallel to the axes (following the ordering of the axes) . As 
n becomes large, the weights on 7r„ will satisfy a strong law of large numbers, 

[n||u2-y||ij 



Let Ui and U2 to converge appropriately on rational vectors to y, use continuity of p(-) 
to finish the proof and verify the limit for all points. 



Continuity at the boundary. Let ek,d = (0,0, ••• ,0,Sk+i,-'' )^d) G R+, and let 
y = (yi, • • • , yfe, 0, • • • ,0) with 7^ if i < A;. Define = y + 6^,^ e int(R^). Then, 




(2.33) 



Symmetric bounds hold for ^[nuij,[nyj- These, along with (2.32) give 



p(ui) < lim n Mog^Lnyj < lim n MogZL„yj <p(u2). 



(2.34) 



p(ye) = hm n ^ log Z^nyej 



d 




d 



p(y)+E(a;(0))/5 ^-S- 



(2.35) 



i=k+l 
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For a reverse inequality, given u e define the following set of paths on 1J^~^: 

Decompose the partition function according to the lattice points where it enters a new 
level in the direction (including an irrelevant weight a;(0) = at the origin): 

^Lnuj = XI n ^(v',i),(v*+i,i)- (2-36) 
For TT e Al„uJ let denote a summand in (2.36). 

l0gZ^= XI log^(vi,i),(vi+i,i) 
1=0 

\nua\ \nud\ 

i=o i=i (2.37) 

i=l 

< log •^|^„u^^^_jj_|_|^„„^e^_^j . 

'^\nud ij+Lnurfed ij ^bove is a partition function in Z*^"^ with weights coupled with the Orig- 
inal weights so that we have the identities 2'(vi,i),(vi+i,i) = ^i+i and /3 <x'(v*, i) = 
loe 

6 [nui,d_ij,[nui,d_ij + [nudjed_r 

Let M — lA^nuj |. Counting the number of ways in which the length from to \nui\ 
can be decomposed into \nud\ + 1 segments gives 

^= TT fL--d + LnMY (2.38) 

By Stirhng's formula, 



M 



{d-i 
n X (i'^i + '^d) log(wi + Ud) - Ui log Ui - Ud log lid^ + o{n) 

= exp {nLu^^_^ {ud) + o(n) } . (2.39) 
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Let £ > and fixed. Assume that — Sd > 0, positive but sufficiently small 
(depending on the value of e). If log Z^nu\ > np{ui ii-i) + n£ there must exist a summand 
in (2.36) with total weight no less than M-ie"(p("i-<*-i)+''). Using this, (2.37) and the 
concentration inequality in Lemma 2.12 [Proposition 3.2.1 (b) in [13]] we compute 

P{ log^Lnuj > np(ui,d_i) + ne] 

< MP{ logZ[;„^^^_^j+L„„,e,_ij > np(ui,d_i) + ne- logM} 

< MP{ logZ[;,^^^j+Ln.,e,_,j > np(ui,d_i + uaea-i) + ns/2 - logM} 

< 2 exp{-n(ci£2 - 2Lu^ (e^)) + o{n)}. (2.40) 

For Ed sufficiently small, Lui^_i(£d) is smaller than e"^. A Borel-Cantelli argument 
gives that 

lim n-^ log Zl^uJ < p(ui,d-i) + £, P - a.s. (2.41) 

n— >-oo 

To finish the proof, consider ye as before. Iterating the arguments from above, 
starting from (2.36), we conclude that 

Pjlog^L^y^j >np{y)+ns} 

< 2'^-^ exp I - n [cis^ - ^ Ly, (e,) + /(e^, • • • , £fe+i)) + o{r. 

where /(e^, • • • , £jfc+i) — > uniformly as the £j — > 0. By Borel-Cantelli we have 

lim n-^ log ZL„y,j < p(y) + £, P - a.s.. (2.42) 
Equations (2.35) and (2.42) give the continuity at the boundary by letting e tend to 



Convergence in all directions. We now show that convergence happens P- a.s. simul- 
taneously for all endpoints y G M^jj. . Note that there is a full probability event of Vl where 
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the convergence is true for all y e Q*^. By (2), it suffices to show that the conclusion is 
true for y with ||y||i = 1. 

Let 5, £ e Q+ and assume < S < e. We define two partitions: 

^tli^^S) = {Vi,V2,...,Vm} 

is a partition of the hyperplane ||y||i = 1 + s, y G M^|_, so that minj^j ||vj — Vj||i < S 
with V e Q*^. Also we project 7r+ to get 

7r]^(£, 5) = {vi - 2£ed, V2 - 2£e<i, . . . , vm - 2£ed}, 

a partition of the hyperplane ||y||i = 1—e. Any point that falls outside the ffist quadrant 
we remove from 7i^{6,6). Observe that for any y with ||y||i = 1, there exist partition 
points Vj e nlj{e, S), Uj = Vj - 2eed e n]^{e, S) with ||vj - y||i + ||ui - y||i < 4£, with 
the vector ordering < y < Vj. For any n e N and choice of vectors Uj, Vj fix a path 
7r„. from [nUjJ to [nvjj . Restrict the space of environments further by assuming the 
following law of large numbers: 

lim n-^ V \uj{xk) \ = ||v. - Ui||E|a;(0)| (2.43) 
Then, following the calculations as in (2.32), 

log ^Lnu^J + log ^Lnu .J , [nyj < log ^[nyj < log ^[nv^J " log ^[nyj , [nv^J (2-44) 

gives 

p{ui) — Ce < lim n"^ logZ^nyj < lim logZ^^yj < p(vj) + Ce. (2-45) 

n— >-oo n— )-co 

Let £ — > along rationals and Uj and Vj to converge appropriately on rational vectors 
to y. Continuity of ]?(■) on rationals gives the conclusion. 
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To prove (c) observe that for any vector y = {i/i, 1/2, ... ,yd) and a permutation on d 
elements a e Sym{d), the following equalities hold: 

p(y) = p(yi, z/2, ■ ■ ■ , Vd) = p(ya{i),ya{2), ■ ■ ■ , ya{d)) = p(o-(y)) 

Restricted on the hyperplane ||y||i = 1, p(y) is concave, symmetric about d~^l. This 
implies (c). □ 

Proof of Proposition 2.4- Without loss of generality, let t — 1. Using the definition (2.5) 
we immediately get the lower bound 

p{d-'l) < (2.46) 

For the upper bound, let £ > 0. Use the same partition 7r^{e) as in the proof of the 
previous proposition and the last part of (2.44), to estimate 



MogZr=n-Mog( Yl l^g^Lnyj) 
y=l|y||i=i 

< log ( n'^ max log Z^ny} ) 
^ y:||y||i=i ^ 



W|v,-y|liJ 

< n~^dlogn + max log ZuvJ + max > (]\u{uj)\, (2.47) 
where Qi = {y ■ ||y — Vj||i < e}. Take n — >■ 00 in (2.47) to conclude that 

p*°*(l) < sup p{\i) + Ce <p{d-H + e,--- ,l + e) + Ce, (2.48) 

||u||l = l+£ 

where C is the limiting last passage percolation constant if the path takes 2ne steps. 
Let e tend to to get the conclusion. □ 

CoroUciry 2.13. There exists a full ¥ -probability event so that 



n 



-1 



log ^[„uJ,LnvJ = P(v - u) (2.49) 



simultaneously for all (u, v) G Ml x Ml with u < v 
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Proof. Prom the previous proposition, the result is true if u = 0. Then we can restrict 
further to all rational vectors u e and use simialr approximations as before to get 
the corollary. □ 

The remaining part of this section is about large deviations results. 

Lemma 2.14. Suppose that for each n, L„ and are independent random variables. 
Assume that the limits 

\{s) = - lim n-MogP{L„ > ns}, (2.50) 
= - lim n-^ logP{Z„ > ns} (2.51) 

exist and are finite for all s e M. Assume that X{a\) — (l){a^) — for some a\, e R. 
Assume also that A is continuous. Then for r e R 



hm logP{L„ + Zn>nr}^< 



— inf {0(r - s) + A(s)}, r>a^ + ax 

0, r <a^^ ax- 

(2.52) 



Proof. The lower bound > follows from 

P{L„ + Z„ > nr} > P{L„ > ns}P{Z„ > n(r - s)}. 

Since an upper bound is obvious, it remains to show the upper bound for the case 
r > + ax- Take a finite partition ax = qo < ■ ■ ■ < qm — '''' ~ 0'4>- Then use a union 
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bound and independence: 



F{Ln + Zn> nr] 



m—1 



+ ^ P{L„ + Zn> nr, nqi < L„ < ngj+J + P{L„ > nqm} 



m—1 



< P{Z„ > n(r - qo)} + ^ ^(^ " <li+^)MLn > nqi} + P{L„ > nqm}. 



Note that A(go) = — ?m) = 0, refine the partition and use the continuity of A. □ 

2.4 Existence of the rate functions 

First a general lemma about existence of hmits of almost super additive sequences. 

Lemma 2.15. Let {onlneN « sequence such that Un+m ^ (in+(im+Cn,m where \cn^m\ < B. 
Then the limit lim^^oo "'t-^^o^ exists (and is potentially infinite). 

Proof. Let 7 = lim^^oo ''^"''^On and assume first that 7 < 00. Specify an £ > and 
let Nq e N such that N^^aNo > 7 — £ and Nq^B < e. We can write any n e N as 
n = kNo + r. Then 



i=0 



From this 



lim n-^ logP{L„ + Z„ > nr} 




n > n ^ttkNo + n — n ^B 



>{n- r)n-^{NQ^aNa - N^^B) + n'^ar - n'^B 



>{n- r)n~^(7 - 2e) + n'^Ur - n'^B. 
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Take lim„_^^ on both sides and let £ — > 0. If 7 = oo, pick any large constant C and let 
A^'o large so that Ng^a^Q > C. Identical steps show that the hmit is infinity. □ 

2.4.1 The fixed endpoint case 

Proof of Theorem 2.5- Existence. For m,n & M+ let x^^^ e {0, l}*^ so that [(m + n)uj = 
[muj + [nuj + „. By super additivity, independence and shift invariance 

P{log^[(m+n)uj > (m + n)r} 

(2.53) 

> P{logZL^uj > mr}P{logZLnuj > nr}P{log Zx^,„ > 0}. 
By assumption (2.9) there is a uniform lower bound P{logZx^ „ > 0} > p > 0. Thus 
t{n) = logP{log Z^nu} ^ is superadditive with a small uniformly bounded correction. 
Similar reasoning shows that either t{n) — —00 for all n or then t{n) > —00 for all 
n> Uq. Consequently by superadditivity the rate function 

JJr) = - lim n-MogP{logZ|„u| > nr} (2.54) 

n— >-oo 

exists for u = {ui, . . . ,Ud) G and r G M. The limit in (2.54) holds also as n ^ 00 
through real values, not just integers. 

Similarly we get convexity of J in (u, r). Let A e (0,1) and assume (u,r) = 
A(ui,ri) + (1-A)(u2,r2). Then 

n-MogP{logZL„uj > nr} > A(An)-MogP{logZL„Auij > nXn} 

+ (1 - A)((l - A)n)-4ogP{logZL,(i_A)u.j > n(l - A)r2} + o(l) 

and letting n — )■ 00 gives 

Mr) < XJuM) + (1 - A)Ju.(r2). (2.55) 
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Similar arguments give existence and convexity of the rate function 

Jt{r) = - lim logP{logZ[°*j > r}. (2.56) 

□ 

2.5 Behavior of the rate functions 

We first need to show that the rate function, whose existence was estabhshed in the 
previous section, is not trivial. We establish nontrivial bounds and the continuity of the 
rate function on the boundary. 

Proof of Proposition 2. 6. Let u = (wi, 1^2, ■ ■ ■ , u^) and let Ui^^ = (wi, 1*2, ■ ■ ■ , tijt, 0, . . . , 0) 
where e {1, 2, . . . , d — 1}, u e R^. We start with an upper bound for J. Fix r e R. 
Superadditivity of the partition functions gives 

{d \nui\ \ 
log^Lnui.fcj + ^{^1,3) ^ \ 
i=k+i 3=1 ) 

d f \nui\ \ 

< Ju,,{r) - V lim n-^ logPi \nui\-^ V a;(a;,j) > \ (2.57) 

' n— >^oo I ' I 

i=fe+l k j=l ) 

d 

= E «^inf^(^)> (2-58) 

i=k+l 

since the sum inside the braces of (2.57) is a sum of i.i.d. uj{xj) r.v. for which the Cramer 
rate function exists by Assumption 2.2. This bound holds as long as 7^ 0. 

For the lower bound, we use the same decomposition as in the continuity of the 
limiting free energy. Recall (2.39) and define for 2 < i < 0? 

i-l 

(u) = ^ (^{uj + Ui) log{uj + Ui)- Uj log Uj - Ui log . (2.59) 
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With this notation, M in (2.39) is given, using a Stirhng approximation by M — 
exp{nFd_i(u) + o(n)} < cxp{nFd-i{u) + ned-i}- The error e^-i can be as small as 
possible for n sufficiently large. The functions Fj_i are jointly continuous in ui, . . . , Ui-i 
and as long as rtj_i (and possibly more of the itj's) tend to 0, so is the function. 

We use this to bound J from below with a standard union bound: 



Juif) > — lim n ^log 7 PjlogZ^r > nr — logM} 

n—^oo ' 

= ^ui,d_i+«derf_i(r - Fd-i(u) - ed_i) - Fd-i{u). 

In the last step above a little correction as in (2.53) replaces [nui_d_ij + [nUdBd-il with 
|_nui,rf_i + nuded-i\. 
Set 

d 

ui,i_i = ui,j_i + E Ukei^i. (2.60) 

k=i 

Proceeding inductively, we get the lower bound 

Ju(r)> Js,,,(r- J2 (^i-i(SM)-ei-i)) - E ^-i("m), (2-61) 

fe+l<i<ci k+l<i<d 

where the e^'s are quantities that come from the errors from repeated use of Stirling's 
formula and go uniformly to 0. Joint continuity of the function F(u) translates to joint 
continuity of the lower bound in (2.61). 

Equations (2.58) and (2.61) suffice to give continuity of the rate functions on the 
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boundary of R^. Let 1 < k < d and assume without loss of generality that the coordi- 
nates that converge to are the last d — k, so that 

with Ui > 0. Then both bounds (2.58) and (2.61) come together and give 

lim JuM(r) = J(„,,...,„,,o,...,o)(r). (2.62) 

When (2.61) is used, it is under the assumption that Ui, . . . ,Uk > 0. Then, the rate 
function Jg^ ^^{r) is continuous in (ui^^, r) by Theorem 2.5 applied when we are restricted 
on the /c- dimensional facet of M!^ and the finiteness assumption of the rate function. 

Assume now that Jo(r) < oo and u^"') 0, u^'") e intM^[.. For an upper bound at 
the origin, one can repeat the steps (2. 57), (2. 58). The partition function now is just a 
sum of i.i.d. a; variables. Then, 

d 

Ju(m){r) < v!r^J,,{r/ut^) + y inf /(x). (2.63) 

1=2 

This, along with (2.61) give that for any (5 > 

lim u]^^ Jei{{r — S)/u^^^) — S < lim J^(m){r) 

< Ju(m){r) < uf'^J^^{r/u]J'^) + 5. (2.64) 



m->oo m— >-oo 



Let ro such that J(ro) — 0. Concavity of Jej gives a further upper bound a;oo(r — wiro) — >■ 
XooT. For a lower bound, take a sequence tn — > x^o- Then, for any 

lii Je,((r - 5)/ur) > tn{r -5)-ui logE(e*"'^W). 

Take lim^^^g and then tn — >■ Xoo- □ 
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Proof of Proposition 2. 7. We prove the second one. Let £ > and N large enough, so 
that En~^ log Z^nu] < + £ for n > A^. Then, 

2exp{-C2(u)£2n} > P{| log^L^uj - Elog^L^uj I > ns} 

> P{log Zl„uJ - E log > ne} 

> P{log Zl„uJ > np{u) + 2ne}. 

This suffices for the result. □ 

Proof of Proposition 2.8. Let m > and n > 0. Let /3 > 0, r G M and let u G M^j. and 
for simplicity assume ||u||i = 1. Observe 

IILnuJIU 

log^f„.j<log(|n(LnuJ)|^max^^( H ^^^^^0) 
< nC + /3r([nuJ) 

Hence 

/3r(LnuJ) < logZf^^j < nC + /3r(LnuJ). (2.65) 
Since is left continuous at r, we can estimate from below 

lim Ju(/3r) = — lim lim n~^logP{logZf , > n(5r} 

> - lim lim n^MogPlnC + ^T(LnuJ) > njir} 

/3->cxD n->-oo 

= lim (r - Cp-^) 

The reverse inequality also follows with the same estimate and the left hand side of 
(2.65). □ 
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2.6 Quenched Large deviations for the polymer path 
and endpoint 

For this section, we restrict on the P - full probability event 

Qo= (cu: lim n-MogZL„uj,L„vj}=Kv-u),V(u,v) e M^[. x R^, u < v| 

whose existence is given by Corollary 2.13. We make no special mention of this fact in 
the proofs that follow. 



Proof of Theorem 2.9. Let ||u||i = 1. We compute 



Z^a;.6n([nul) ll7=l ^ 



-n-HogQ^^{xn = [nu]} = -n'Mog -^--^"([""1) A^^' 



^xo,„entot(n)lij=i'^ 

—n~^ log 



Z[nu] 



Ztot 

-p(u). (2.66) 

□ 

Proof of Theorem 2.10. We only treat the free endpoint case. The constrained endpoint 
result follows by similar arguments. 

Fix L e R+ and let 7 : [0, 1] — >■ be a curve such that each coordinate 7j(i), 
1 < J < c^, is non-decreasing and L-Lipschitz. Since 7 is Lipschitz, it has a derivative 
almost everywhere. For the upper bound, this is the only fact of 7 that we are going to 
use. 

Pick £ > and jVeIt) an £— neighborhood of 7 in the || • ||i norm. (For the definition of 
this neighborhood, consider 7 as a set in R^.) For this choice of £ > 0, let M sufficiently 
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large to define a partition of the time interval [0, 1] 

ttm = {0 = < ^2 < • • • < = 1 : 7j(^m) - ^j{ti) < £/2rf, 1 < i < d}, 

and assume without loss that for all partition points ti, ^'{ti) exists. 
Abbreviate 7(^1) = Uj and define the rectangles 

R{\nUi\, \ nni+i\) = {v e Z"* : [nujj < v < [nUj+iJ, 

^ ||v - [nujj lloo < II \n\Xi+i\ - [nuj ||oo}. 

The definition of ttm implies that the disjoint rectangles i?([ujj, [u^+ij) C nMe{'^)- 
Microscopically, for any path a;o,n define by Xq ^ to be the piece of the original path that 
lies in i?([nUjJ, [nUj+iJ). To get an estimate for the quenched rate function from above, 

Qn{^i),n e nNeil)} > Qn{xo,n gocs through [nuij , [nuaj, • • • L^UmJ.} 

= Qn{^o,n ^ L«Ui+iJ),4,„ Starts at [nuj, 1 < i < M} 

nM-l y 
_ i=l ■^[nuij,[nui+ij 



Take logarithms in (2.67), divide by —n and let n — > 00 to conclude that 



(2.67) 



M-1 



lim n-' log Q'^{xo,n e nKil)} < p(d-'l) - Pili^i+i) - li^i)) 

n— >-oo . 1 

1=1 

M-1 „f.^^ 

^p{d-H)-J2p{ ii^)ds) 

i=l 

M-1 „t,+i 

<p{d-'l)-J2 / P{l'{s))ds, 

i=l 



(2.68) 



where (2.68) is the result of Jensen's inequality applied on the concave function p and of 
the fact that p is 1-homogeneous, by Proposition (2.3). Equation (2.68) gives the upper 
bound. 
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Figure 2: The idea of the proof of the upper bound in Theorem 2.10. The rectangles are 
defined using the the partition points (the circles). They are disjoint and inside the e— 
neighborhood, so any polymer chain inside those rectangles is inside the neighborhood. 




Figure 3: The idea of the proof of the lower bound in Theorem 2.10. The diagonal lines 
are the ||7(ti)|| hyperplanes. The extra error (compared to proof of the lower bound) 
comes from the common small rectangles with the hyperplanes as diagonals. That error 
can become small if the e-neighborhood is sufficiently narrow. 
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The remaining proof is about the lower bound. A word on the arrangement of the 
proof: The first part of what follows is used to specify an > such that certain 
macroscopic uniformity and continuity conditions are satisfied (here we use the fact that 
7 is Lipschitz). After Eq is specified we define a neighborhood for e < Eq and work 
microscopically to bound the quenched probabilities. That part of the proof works for 
any e < eo- The need for the conditions that specify Eq become apparent after taking 
the limits (calculations (2.75)-(2.77)). 

Fix L e M+ and let 7 : [0, 1] M!^ be a curve such that each coordinate 7j(i), 
1 < J < c?, is non-decreasing and L-Lipschitz. An immediate consequence is that 
< Jj{t) < L for all t where the derivative is defined. Since 7 is Lipschitz, the difference 
quotients are bounded 

^^^^^^^|<L x,ye [0,1], 1 <j <ci. (2.69) 

Let s' > and restrict the limiting point-to-point free energy p(-) on the set ^4^ = 
{||u||i < dL}. The function p is uniformly continuous on ^4^, so it admits a modulus of 
continuity and we can specify S' > so that 

sup \p{u) - p{v)\ < e' . (2.70) 

||u-v||i<5' 

For a given s' > 0, denote by n(e') the supremum of 5' > such that (2.70) holds. As 
soon as S' is specified, define a partition 

Am(£') = {0 = ti < t2 < • • • < = 1} 

so that for all partition points ti, the derivative j'{ti) exists, and so that for all 1 < i < M 
and 1 < j < d, 
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Pick £ > and N'eil) a-n £— neighborhood of 7 in the || • ||i norm. Assume e is small 
enough so that 

6 <eo = mm{Q{6'/2M), 6'/2M, min {\\-f{U+i) - -f{ti)\\i/M}} . (2.72) 

L l<i<M— 1 J 

For any vector u G M.'^ define the hyperplane H{u) = {x G M^j. : ||x||i = ||u||i} and 
the positive half-space -f^+(u) = {x G : ||x||i > ||u||i}. We also denote -f^-(u) = 

Observe that there exist vectors {a^}i<i<M-i such that for each i,j the following 
hold: 

1. B.} < -i{U) < af , l<i< M, 

2. ||7(i,)-a^||i = £, i = l,2, 

3. R{[n3.]l Lriaf+iJ) D nMMUA+i)) n H+{[n^{U)\) n H_{[n^{U^,)\) = . 

The set l-L^ is the set between two consecutive hyperplanes and in nMe (see Figure 
3) . Denote the partition function by 

[n7(ti+i)J 

z«"= E n 

where x^*^ is any up- right path living in "H". 

Define = {xo^ne ■ = L^a^^JI^ and nf„^ = {xo^ne ■ a^o e H{[nj{ti)\),Xne = 
[na^J}, the set of all path of length ne that start from [najj and the set of paths 
that start somewhere on the hyperplane H{[n'y{ti)\) and end after [ne\ + 0(1) steps at 
|_na^J , respectively. The error comes from the integer parts, but is eventually immaterial 
so we ignore it for convenience. 
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Observe that any polymer chain that starts at [na^J and ends at [riaf,,.!] has to cross 
the hyperplanes H{[n'y{ti)\) and iJ( [n7(tj+i)J ) at some points a and P respectively. 

The three conditions, along with (2.72) guarantee that Ui<i<M-i L^^lJ > L^^i+iJ ) 2 
nMe{'^) while the common i?([naJJ, [nafj) (the red-shaded rectangles in Figure 3) are 
pairwise disjoint and "small", in the sense of the following bound: 



(a,/3)ei? ( Ln7(ti)J ) ( Ln7(t,+i)J ) 

[nej [nej 



\ne\ [nej 

>ZHn min JJe'^'"^^^-) min JJe^'^^^^^. (2.73) 
Hence we can bound the probabilities 

a;o,neAre(7) j=l 
M-1 

< (z„)-i J] 

i=l 

M 

lne\ [nej _^ 

xf min TTe^'^^^^^x min TT e'^'^^^^)) . (2.74) 

We are going to take logarithms on both sides of (2.74). On the right-hand side, we have 
logmin^-jj^^gni ^ 11^=^1 e~^l'^*^^J^I which is the first passage percolation time of ne steps if 
the weights were negative (so the negative of last passage percolation time if the weights 
where positive) with limiting constant c_\i^\{e). 
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Then, we conclude 

M-l 

- EE logQ:^{xo,„ e nA4(7)} > - J] p(a?+i - a,^) + 2Mc_|,|(£) 

i=l 

M-l 

i=l 
M-l 

> p{d-'l) - J2 Pbiti+i) - liU)) - Ce' (2.75) 

i=l 

Pid V-Z^P[ — t+i-t- — jiti+i-ti) -Ce 

(2.76) 

> p{d-H) - ^ piiimu+i - U) - Ce'. (2.77) 

i=l 

Equation (2.75) is the result of (2.72) and (2.70). Homogeneity of p gives (2.76) and 
then (2.72) and (2.71) give (2.77). The bound in (2.77) is independent of e so letting 
£ — > does not affect it. To get the result, let the mesh of the partition tend to and 
then let s' ^ 0. □ 



M-l 
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Chapter 3 

The log-gamma Polymer Model 

3.1 The model 

3.1.1 Introduction and results 

In the present chapter, we derive an exphcit expression for the upper tail rate function 
for log Z^nu\ in the case of the 1 + 1 dimensional log-gamma polymer. The computations 
are tractable exactly because of the Burke property More detailed information about 
the model and basic properties can be found in Section 3.1.2. 

The results for the particular 1 + 1 dimensional log-gamma model. The distributions 
of the uj weights are i.i.d. 

~ logy, where ~ Gamma(/x), (3.1) 

where the density of the Gamma(/i) is given by (A. 2). For this choice of i.i.d. weights, 
denote by 

M,{i) = logE(e«-(°'°)) = logr(/. - - logr(/x) (3.2) 

the logarithmic moment generating function. 

It is convenient for the proof of these results to write the vectors in using both 
their coordinates. The main result is an explicit formula for the upper tail large deviation 
rate function for the logarithm of the partition function. 
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Theorem 3.1. Let r e R and let uj{i,j) to be distributed as in (3.1). The function 

Js,t{'^') defined by (2.13) is given by 

Js,t{r)= sup -inf (tMeii) - sM^.e{-^))] ■ (3.3) 

Remark 3.2. While the symmetry of Jg^t 'is clear by definition (2.13), it is not imme- 
diately obvious from (3.3). From the proof of the theorem for s < t one can check that 
we can restrict the set where the inner infimum is taken to 9 & [{/i + ^)/2,//). Under 
the assumption that s > t it is easy to check that the infimum is now attainable for 
9 e (^, (// + ^)/2]. If s < t and 9q gives the infimum 7*, it is easy to check that in- 
terchanging the roles of s,t will give the same infimum ^y^ at 9i — /i -\- ^ — 9o. These 
symmetries will be explained in detail in the proof the theorem. 

The next result is about the free-endpoint case. 

Theorem 3.3. Let s > 0, r e M. The large deviation rate function (2.14), for /3 — 1 is 

- lim n-^ logP{log > nr} = Js/2,s/2{r) (3.4) 

n— >-oo 

where Js,t{f^') is given by (3.3). 

3.1.2 The model and the Burke property 

For the rest of this chapter we also the the parameter /3 to equal 1. Henceforth we adopt 
the multiplicative notation for the polymer measure. It is also convenient to adjust the 
notation for the partition functions and redefine the rate functions: 
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On each site of we assign weights Yij. For any {k,£) < (m,n) define the 
partition function for paths that start from {k,i) and whose endpoint is constrained to 
be {rrijn), by 

m+n 

where Il{k,i),{m,n) is the collection of up-right paths x. — {xj)k+e<j<m+n inside the rect- 
angle -R™/' = {k,k + l,...,m} X {i,...,n} that go from {k,i) to (m, n): Xk+e = {k,i), 
^m+n — {rn,n). We adopt the convention that ^(jt,^),(m,n) does not include the weight 
at the starting point. In the case that the weight at the starting point needs to be 
considered we also define 

m+n 

Zfk,e),{m,n) = X] n -^^i ^ yk,iZ(k,e),im,n), (3.6) 

If a value is needed, then assume that Zk/^{k/) = 1. In the special case where (/c, t) = (0, 0) 
we omit the subscript from the above notation and we also set Fo,o = 1- 

We assign distinct weight distributions on the boundaries N x {0}, {0} x N and in 
the bulk N^. To emphasize this, the symbols U and V will denote the weights on the 
horizontal and vertical boundary respectively: 

Ui,o = Yi^o and Vqj = Yqj. (3.7) 

Our results depend on the explicit distribution of the weights; all weights arc recip- 
rocals of gamma variables. To be precise, here are the assumptions on the distributions: 

Assumption 3.4. Let < 9 < < oo. We assume that the weights {Uifi,VQj,Yij : 
i,j&N} are independent with distributions 

U^Q ~ Gamma{9), Iq~-^ ~ Gamma{fj, — 9), and Y-^^ ~ Gamma{fj,). (3.8) 
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Given the initial weights {Uifi,Voj,Yij : i,j e N} and starting from the lower left 
corner of N^, define inductively for (i, j) e 

The partition function for the model with the boundary condition is denoted by Zm,n 
satisfies 

= YmAzt\n + for (m, n) e (3.10) 

and one can check inductively that 

jj _ ^m,n A ^/ — "''^ 1 1 ^ 

^m,n — — ana Vm,n — " (^) • W-llj 

m— l,n m,?i— 1 

Equations (3.10) and (3.11) are also valid for Z^{m,n) since the weight at the origin is 
canceled from the equations. 

The key result that allows explicit calculations for this model is the Burke-type 
Theorem 3.3 in [29]. 

Let z. — {zk)kez be a nearest-neighbor down-right path in Z^, that is, Zk G and 
Zk — Zk-i — ei or —62. Denote the undirected edges of the path by fk = {zk-i, z^} and 
let 

Uz^ if fk is a horizontal edge 

if fk is a vertical edge. 
Let the lower left interior of the path be the vertex set X = {(i, j) G : 3 m so that (i-|- 
mj + m) e {zk}}. 

Recall the definition of Xij from (3.9). 

Theorem 3.5 ([29]-Burke Property). Under the assumption (3.4), and for any down- 
right path {zk)kez in 2+? the variables {Tf^,Xz : k E Z, z E 1} are mutually independent 
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with marginal distributions 

~ Gamma{9),V~^ ~ Gamma{ii — 9), and ~ Gamma{ii). (3-12) 

Prom this, one can compute 

E(log Zj^l) = mE(log U) + nE(log V) = -m^o{d) - n^o{l^ - 0). (3.13) 

In [29] a law of large numbers is proved for the limiting free energy in the case of no 
boundary weights. Let (s, t) e and observe that there exists a unique 9 = Og^t £ (0, A*) 
such that t-^iin -9)^ Define 

p^{s,t) = -{s^o{9s,t) + i*o(/^ - 9s^t)). (3.14) 

For the model without boundaries with Y~^^ ~ Gamma(/i), i > 0, j > 0, the limiting 

free energies can be evaluated explicitly to be 

lim n"MogZL„5j Lntj =Pm(s,^) P - a.s. (3.15) 

n->oo 

and 

lim n-^ log^|°* I = -s*o(a*/2) P - a.s. (3.16) 

3.2 Continuity of the rate function on the boundary 

The logarithmic moment generating function of the bulk weights log j ~ — log Gamma(/x) 
is 

iogr(/i-0-iogr(//), e</^ 



M^(0 = logE(e«i°^^) 



oo, C > /W- 
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Its convex dual is the Cramer rate function (recall (2.8)) 1^ of this distribution, given 

by 

/^(r) = -r^-i(-r) -logr(^o^(-r)) +/ir + logr(/x), r G M. (3.17) 

From Theorem 2.5 we know that the rate function for the model with i.i.d weights is 
continuous on the boundary, and in this case, its value equals the Cramer rate function 
for the sum of i.i.d. weights with inverse Gamma(/i.) distribution. 

The rate function on the boundary is given by 

xlf^{rx~^), r > — x^o(/^), x > 

Jo,xi^) = Jx,o{^) = ^ r/x, r > 0, X = (3-18) 

0, r < -x*o(a«), x>0. 
For the second branch of (3.18) we used 

limxlogr(*o^(-rx-^)) = 0, r > 0. (3.19) 

Note that Jo,o{f) — r/J, as given by (2.15). (Here // = x^o)- One can also interpret the 
second branch as the large deviation rate function for a single logF random variable, 
where ~ Gamma(//). 

The strong law of large numbers for the limiting constant of the free energy from 
(3.15) and continuity of J give the support described in the theorem. The fact that 
J is strictly increasing for r > PfM{s,t) can follow independently after finding the rate 
function explicitly in Theorem 3.1. 



3.3 Exact point-to-point rate function 
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Lemma 3.6 (Varadhan's lemma). Let Z^nsi,[nti the partition function given by (3.5) with 
weights Y such that ~ Gamma{ii) . Also let Js,t{f) the upper tail large deviation for 
log Zi 

ns] , [ntj • 

Then /or < ^ < 

lim n-MogEe«^°s^L-j,Lntj = sup{r^ - Js,t(r)} = J*tiO- 



Proof. Let 7inf = lim„^^ n ^ logEe^^°^^L'»sJ.L"tJ and 7s„p = lim„^oo n ^ logEe^^°^^L"sJ,L"tJ 
First an exponential Chebychev argument for a lower bound: 

n-HogF{\ogZlr^s},lnt} > nr} < -^r + n~MogEe^'°s^L-J'L"*J. 

By letting n — >■ oo on a suitable subsequence we get that for all r e R 

7inf >^r- Js,t{r)- 

Optimizing over r we get 7inf > Js,t{0- 

For the upper bound, first note that there exists a > 1 such that < jji, 

sup [Ee^^^^s^L-J'LntjV^" < oo. (3.20) 
To see this, distinguish cases where < 1 or otherwise. For cu^ < 1, 

l/n Lf^sJ + HJ ^ * - - 

xsn.([ns\,[nt\) j=l 
[ntJ + \ns\ 

<{ n 



a;en([nsj,[ntj) i=l 
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For > 1, we use Jensen's inequality. Let N denote the number of paths. 

l/n L"«J + L"*J 



a;en([nsj,[ntj) j=l 



l/n 

a;en([nsj,[ntj) i=l 



Now, we can show that 

hm EE n~MogEfe^^°s^L-j,Ln«ji{iogZ|„,| > nr}) = -oo. (3.21) 
Using Holder's inequality, 

-MogE(e«^°^V.J.L"*Jl{logZL„,j,Lntj > nr}' 



< Q;-Mogsup(Ee"^'°s^L-J'L"*J)V"+ ^^^^-ilogP{logZ|„,| i^ti > nr). 

Taking a limit n — > oo we conclude 

^ n-MogEfe«i°^^L-j,Lntji|iogZL„,j,L„tj > nr}) < C - ^^J,,t(r). (3.22) 

Letting r to oo concludes the proof, since JsA"^') ^ non-constant increasing convex 
function. 

To estabhsh an upper bound let 5 > and partition R so that for i e Z, r, = ib. 
Then for any m 

oo 

Mog (^Ee^i°g^L-j,L-j j = n^Mog ( ^ Ee«'°s^L-J.L"*J l{nn < log Z^^^j.^ntj < ^n+i}) 

i=— oo 
m 



n 



<n-Mog( ^ e^^'-^+iPjlogZL^.j^Lnij > nr,} 

+ E(e«^°s^L-J.L-Jl{logZL„,j,Lntj > nr^})) 



(3.23) 
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A limit along a suitable subsequence yields 



a — 1 



7sup < max { max {Cr^+i - Js,t{^i)}^ i^-m, C Js,t{rm)}] 

—m<i<m (X 

a — 1 

< max{sup{^r - Js,t{r)} - S^, Cr_^, C JsA^m)}} 

To finish the proof, let 5 — > and m — > oo. 



□ 



An immediate consequence is that logZ'-' and logZ have the same rate functions. 
This follows from having the same convex dual for < ^ < (for C > A* both rate 
functions are oo): 

lim n~MogEe^^°^^°"''J.Ln*J = lim n~MogEe^^°s^Ln'.J.L"tJ+«i°g^i.i 

n->oo n->oo 

= j;,(0+ lim n-HogWYl 

n— >oo 

= KM)- 

We are using this fact in the remaining part of the section without alerting the reader. 

For what follows we need some notational conventions. Assume the polymer lives in 
the rectangle {0, . . . , \ ns\ } x {0, . . . , \ nt\ } and let — \nt\ < k < [ns\ . Let Yq^ ~ 1 and 
define 



( lnt\ 

H V,-\ for - [nt\ <k<-l, 

j=-k+l 



V-1, 



(3.24) 



Vo Y\_ ^»,o, for < /c < [nsj , 

1=1 



and 



-[nt\ <k<-l, 



k=<j(l,l), ^ = 0, 

{k, 1), < A; < [ns\, 



(1, [-naj), -t<a<0, 



and M = <{ (1^1)^ a = Q, 

{[na\,l), 0<a<s. 



(3.25) 



a = lim n ^ [naj = < 

n— >oo 
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It is going to be notationally convenient to assume that k — [na\ for some a e [—t, s]. 
Whenever this happens, we identify k = [naj and we assume that n is large enough so 
that [naJ ^ 0. When we take the hmit as n ^ oo to compute the various rate functions, 
we will need a continuous (and scaled) version of (3.25). For this reason we abuse this 
notation by writing 

(0, -a), -t<a<0, 

(0,0), a = (3-26) 
(a, 0), < a < s. 

Observe that for all a e [— the r.v. logT^L^aj is a sum of independent log-gamma 
random variables: For a < 0, log rjina} is just a sum of i.i.d. log Gamma(//— ^), the Cramer 
rate function exists and so are the lower and upper large deviations rate functions. 
For a > 0, 

[nt] [naJ 

log Vina} = ^oT/ - ^i^O- 

0=1 i=l 

Setting Ln = Y^^j^i logK)7/' ^ ~ 1^1=1^ log[/j~Q^, we appeal to Lemma 2.14 so the 
upper large deviation rate function 

- lim n"MogP{log77|„aj > nr) = ^^(r) (3.27) 

n->oo 

exists and is convex and continuous. With these definitions we can have the following 
computational lemma that we are using throughout. 
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Lemma 3.7. Fix a e [—t,s\ and let Ka{r) defined by (3.27). Then 

(i + a)(logr(//-^ + 0-logr(//-^)), -t<a<Q, ^ > 0, 

t ( log r(/. - e + - log r{i,-e)) + a{ log r(^ - o - log r(^)) 

0<a<s, 0<^<^, 

oo, otherwise. 
Proof. Fix a e [— s]. For — i < a < 0, the first branch of (3.28) is the logarithmic 
moment generating function for log-gamma weights. The second branch comes by taking 
the limiting logarithmic moment generating function for ^ > of rj^na] ■ D 

Recall (3.25). Let (a, 6) G [— s]^ and let m^^^a and mj^b the zeros of Ha and J(^s^t)-h 
respectively. Define 

Hs;tii-) = Jhn n-MogP{log77L„aj +log^P„bj,(Ln.j,Lntj) > rir} 

0, r < m^n + mjh 



(3.29) 



inf {tia{x) + J(s,t)-b(r - x)}, otherwise. 

^"i«,o<a;<r-mj,6 



where Ka{x) is given by (3.27). The existence of Hg'f{r) is estabhshed by Lemma 2.14 
and continuity in the b argument (when a, s, t, r are fixed) follows directly from the 
continuity of J in 6, the fact that J is always finite and that x can be restricted in a 
compact set. In the case where a — b we define H^j-{r) = (r). 

Let s,t > 0, r & M.. Let {^j,[ntj }i<j<|^„5j be the weights as defined by (3.9). Define 



( L'' 



R,{r) = - lim n-HogF{ Vlogt/^.Lntj > 

I i=l 



nr 



sleirs-^), r > -s*o(^), 
0, otherwise, 



(3.30) 
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lnt\ 




1 



k 



[ns\ 



Figure 4: Two possible polymer paths with weights described by (3.8). The shaded 
parts explain the decomposition of the partition function Z^nsi,[nti needed for Lemma 
3.8. 



with convex dual 



slogT{9-0-slogT{9), 0<^<9, 



oo, otherwise, 



(3.3i; 



where we use (A. 4) and (3.17) to obtain the first branch. 

Lemma 3.8. Let s,t > , a G [— t, s] and r G M. Let (for fixed s,t,a) Ha, H^^ and Rs 
be defined by (3.27), (3.29) and (3.30). Then 



R,{r)= ini {Hl,ir)}= inf inf {^^(x) + J^,,)-.ir - x)}. (3.32) 

—t<a<.s — t<a<s rriKjaS^Sr—mj a 



Proof. We start by decomposing Z^^^^ \nt\ according to the exit point of the polymer 
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path from the boundary: 

[nsj + [nt] 

^■^^(0,0),([ne],lnt]) 

[nt\ ( / ^ \ ■) L^^'J r / ^ \ "1 

^ U n ^O'^' ) ^(°i/),(M,LntJ) f + U n ^^'0 ) ^{k,lUlnsi,lnti} \- (3-33) 

e=i ^ ^ j=i ^ ^ k=i ^ ^ i=i ^ J 

Dividing both sides of (3.33) by Z^^^^^ — Hj"*! ^oj and through multiple uses of 
(3.11), we get 

[nsJ [nt] , , Ynt\ 



n = ] ( n ^o7/ )^a,^),(Ln.J,LntJ) \ 

1=1 e=i ^ ^j=e+i ^ ^ 



Yns\ , . [ntj k X 

+ 1] 1 ( n ^ot/ n ^^'0 ) ^8,i),(M,Lntj) 

k=\ ^ ^j=l i=l ^ 

[nsJ 

= X/ ^'^^MLnsJ.LntJ)- (3.34) 
k=- [ntJ 

Since all terms are nonnegative we get the bounds 

[nsJ [nsJ 
«=1 fe=-[ntj 

<Ln(s + t)J max r7fcZg(L„,j_LntJ)- 

Take logs of these inequalities to get 

\ns\ . \na\ . 

log r]k + logZ£(LnsJ,LntJ) < X] Vi^Ynt\ < log <^ XI ^fc^k,(Ln.J,LntJ) f 

- , .^^^1 , i + ^k,(Ln.j,Lntj)} + log(n(s + t)). 

(3.35) 
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For any a e [—t, s], 

V'logC/j,L„tj > nr \ 

A — / '1- J 

1=1 ^ 

- Jm^^~^logP{log77L„aj +logZp„^j_(L„,j,LntJ) ^ 

> -HI,{t). (3.36) 

Equation (3.36) is valid for all a, so optimizing over a e [—t, s] gives 

{[nsJ . 
^i.L"*J > "''^ r ^ - inf inf {na{x) + ^(5,t)-a('^ - x)}. 
' I — t<a<s aK<x<r—aj 

i=l 

This proves the upper bound in (3.32). 

The remaining of the proof is about the lower bound. Let £ > 0. Fix a sufficiently 
small 5 > and let — i = Oq < Oi < • • • < = < • • • < = s be a partition of the 
interval — [t, s\ so that |aj+i — a^l < (5. For a given ^ > assume 5 = 5(£, ^) is sufficiently 
small so that b /^"^^ ^0(2;) < e/4. 

Without loss of generality assume Oj > 0. For any integer k e [[riOiJ, [noj+ij], we 
can estimate 
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P{l0g77fe + l0gZg(L„.J,LntJ) >^^} 

log77Lna,+iJ+log^P„ad,(Ln.J,Lntj)- E ^^g C/i,o - ^ loglS-i>nr| (3.37) 

i=k+l j=[nai} 

< P{ l0g77Lna,+iJ + l0gZp„^.j,(L„,j,LntJ) > ^(^ - ^)} 

{[nai+i] k-1 >| 

i=A;+l j=\nai} 

< P{ logr^Lna^+iJ + l0gZp„^^j_(L„,j^LntJ) > ^(^ " ^)} 

< P{ logr?Lna.,.j + log^&,j,(L..j,Lntj) > ^(^ " ^)} + e-l^^-^^/^ *°(^) '^^) 

< P{ log77Lna.+ij + log^P„a.j,(M,LntJ) > ^(r - s)} + e-"«^/^ (3.38) 

The key to the bound is the fact that the upper tail large deviation rate function for 
sums of logF^^ random variables has unbounded slope. This is not true for logF. This 
is why Tjk changes to r]ina,+,i and ^£(L„,j,Lntj) ^°naij,(M,LntJ) (3-37). For the case 
< 0, the corresponding changes will be rji. to rj^nai} because the weights are reciprocals 

and ^k,(M,LntJ) ^°na,j,(Ln.j,Lntj) before. 

Now for the actual error estimate. Assume n is large enough so that ne > log{ns+nt). 
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Equation (3.35) implies 

[nsj 



n 



[ntj > nry 
^ i=i ^ 



[nsj 

<n'Hog P{ log r]k + log ^k,( L„sj , [ntj ) > ^'^ - log(ns + nt) } 

fe=-[ntj 

< max max n^MogPj logr^^ + logZwi^sl l„tn > n(r - e)} 

+ log(ns + nt) 

- o<K^-i''"' (^"f log77L„a,+ij + logZ°^j_(L„,j,L„tj) > n(r - 2s)} + e'^^^/^^ + e 

(3.39) 

- oJ<^m-i''~' ^""^^^ log77Ln„,+ij + log Z° ,^j,(L„,j,L„tj) > n{r - 2s)} V (-^6/2) + 2s. 

(3.40) 

Equation (3.39) follows from (3.38). Take a limit n — > oo in equation (3.40) to conclude 
-Rs(r) < max {-//^r '"'(^ " 2^)} V + 2s (3.41) 

0<i<m— 1 

< max {-<r^''^'(r - 2s)} + 2s (3.42) 

0<i<m— 1 

< max {-H^\+"'''+\r -2e)} + e' + 2s (3.43) 

0<i<m— 1 ' 

< max {-<r'«'+^(r)} + s' (3.44) 

0<i<m— 1 

< sup {-//«^(r)} + s'. (3.45) 

—t<a<s 

Equation (3.42) is the result of ^ tending to infinity in (3.41) and noting that ff^'+i'"' < 



oo. 



Equation (3.43) requires explanation. Observe that for r, s, t, s fixed there exists 
a compact set K that depends only the fixed parameters, that contains all intervals 
Ka,b = [iTT'K,a,^ — 2s — mj^h], for which the x variable ranges over in definition (3.29) for 
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H'^f{r — 2s). This follows from the continuity of m^^a and mj^^ in a, h respectively, and 
from the fact that a, h range over compact sets. Note that by enlarging the compacts 
Ka,b to K we do not change the value of H'^'^{r — 2e). 

For X restricted in K, J(s j)_b(r — 2e + x) is uniformly continuous in (6, x). Then, for 
£' > we can assume the mesh of the partition is small enough so that for any fixed x 

Hence, in (3.43), e' is the error that comes from changing the superscript in H from 
to Cj+i. Equation (3.44) comes from letting £ — > and the continuity of H in the 
r variable that follows from arguments similar to those that justified (3.43). Letting 
£' — >■ in (3.45) gives the lemma. □ 

We will also need the following lemma. 

Lemma 3.9. For a fixed ^ e [0, //) the function 

|-4,t)-a(0, forQ<a<t 
G,{a) = (3.46) 

I oo, otherwise 

is convex, lower semi- continuous on R and continuous on [0, i]. In particular, G*^{a) = 
G^{a) for a e M. 
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Proof. To show convexity on [0, t] let A e (0, 1) and a — Aoi + (1 — A)a2: 

-J*t,t)-^{0 = -SUp{Cr - J(t,t)_a(r)} 

< inf inf {A(J,_„,,(rO - ^ri) + (1 - XM-a,,^) - ^r,)} 

r-eIR (ri,r2): 

Ari+(1— A)r2 =r 

= inf {A(J,_„,,(ri) - eri) + (1 - A)( J,_„,,(r2) - ^r^)} 

(r-i,r2)eIR2 

= A inf { J,_„,,(rO - J + (1 - A) inf { J,_„,,(r2) - ^r,} 

= -XJUAO - (1 - A) (3.47) 

The inequahty comes from the convexity of J in the variable {t — a, t, r). 

For finiteness on [0, t] it is now enough to show that G^{a) is finite at the endpoints. 
Continuity then follows in the interior (0, t). First take a — t. Then Jq ^ is the dual of a 
Cramer rate function, and for ^ > 

Gdt) = -Jo,tiO = -UogEe«'°^^^-° (3.48) 

which is finite for ^ < 

Convexity of Js,t{r) and symmetry Js,t{r) = Jt,s{r) imply Jt,t{r) < Jo,2t{r). From 
this 

Q(0) = -j;,(O = infM,t(r)-M 

< inf { Jo,2t(r) - ir} = - Jo*,2t(C) < oo. (3.49) 

Continuity at a — and a — t. Case 1: a = 0. To show that is also continuous 
at the endpoints, we first obtain a lower bound. For any r e M, 

JUtiO >r^- Jt-a,t{r) 
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hence, by continuity of Jg^t in the (s, t) argument, 

nmj;_„,(0>re-JM(r). (3.50) 

Optimize over r to conclude hma_^o Jt-a,tiO > Jt,t(0- 

For the upper bound, let n e N. Then we estimate, using Lemma 3.6 

= lim n-MogEe«'°s^L"*J.Lntj 

' n— >-oo 

> lim n"MogEe^'°s^L"(t-a)J.L"*J + lim n"MogEe^^*=* "(*-«) J +i^°^^*'L"tJ 

n— >-oo n— >-oo 

= ^t-a,t(0 + «logEy«. (3.51) 

Taking a — > yields the result. 

Case 2: a = The lower bound follows as in the previous case. For the upper bound 
we use a path counting argument. Consider first the case where < ^ < 1. Then, 

[nt\ + [n{t-a)\ ^ 

^r-a,t(0 = linin-MogE( H ^) 

xeU{ln{t-a)i,[nti) i=l 
[nt] + [n{t-a)} 

< lim n-Hog y TT E(r)« 

x€n{[n{t-a)},[nt}) i=l 

= F{t-a,t) + {2-a/t)J*^,{0 
For 1 < ^ < /X, Jensen's inequality yields 

JUAO < m -a,t) + {2- a/t)J*,{0. 

Let a ^ t to get the result. 

The fact G^* — is now a direct consequence of lower semicontinuity, by[24, 
Thm. 12.2] □ 
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Proof of Theorem 3.1. We begin by expressing the explicitly known dual from 
(3.31) in terms of the unknown function J(^s,t)-&- Recall that a e [— t, s] is the macroscopic 
exit point of the polymer chain from the boundary and a is given by (3.26). 

Fix < < ^. By (3.29) and (A.5) we can write Hl^{r) = («;„□ J(s,t)_a)(r). Then 
by (A.4) and (3.32) 

Kii) = sup sup [ri - {KaUJ^s,t)-a){r)] 

—t<a<s r 

= sup (fi;„nj(^,t)_a)*(C) 
—t<a<s 

= sup {<(0 + ^(W0} by(A.6). (3.52) 

—t<a<s 

Equations (3.31) and (3.52) give 

s\ogT{e-i)-s\ogT{e)^ sup K(0 + ^(Va(0}, o<e<^. (3.53) 

—t<a<s 

Define 



u{e) = < 



-h^{e) = M^_,(-o = iogr(/x - e + - iogr(/x - ^), < a < 



d^{e) = M,(0 = iogr(^ - - logr(^), < a < s 

(3.54) 

and substitute (3.28), (3.54) into equation (3.52) to get 

We now speciahze to the case s = t and we treat ^ as a variable in the remaining 
part of the proof. We assume that e (^, //) and ^ > is fixed. When s — t, symmetry 
of Js,t allows us to write (3.55) as 

t{d^{e) + h^{e))= sup {amax{/ie(^),de(^)} + j;_„,(0}. (3.56) 

0<a<t 
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Assume first that (//+^)/2 < 9 < fj,. This is equivalent to h^{0) > d^{0) and equation 
(3.56) turns into 

t{d^{e) + h^{e)) = sup {ah^iO) + JUtiO}. (3.57) 

0<a<t 

Let h^{9) = V and G^{a) = —Jt-a,tiO- This notation makes (3.57) 

t{{d^ o hj^) (v) + v) = sup {av -G^{a)} = GUv), /i^(^) < v < +oo. (3.58) 

0<a<t 

Now assume that ^ < 9 < (/x + ^)/2. Then, equation (3.56) becomes 

t{d^{e) + hm - sup {^d^{0) + J;_„,(e)}. (3.59) 

0<a<t 

LetV'^,^ : (^, (//+^)/2] — )■ [(/x+^)/2, /x), 9 ^ — A*~^+C is a homeomorphism 

between the intervals (^, {ji + ^)/2] and [(// + ^)/2,^). It has the following properties: 
First, d^{9) = h^{l^-9 + () = h^{ipf,^^{9)) and d^{n-9 + ^) = d^{ip^^^{9)) = h^{9), hence 
it fixes the sum 

h^{9) + d^{9) = h^{^,,^{9)) + d^{ij,,^{9)). 
Equation (3.59) can then be re-written as 

t{ddi^,,m + = sup {ah^{^,M) + j;_„,(e)}, (3.60) 

0<a<t 

where h^{'ipfj,,^{9)) G [h^{^^), oo). This shows that equations (3.59) and (3.58) are equiv- 
alent, so we can restrict 9 e [^^,A*) without loss of generality. We will work only with 
equation (3.58) from now on. 
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We compute 

J{t,t)-a{r) = sup {rC - J^aAO} by (A.7) and Theorem 2.5, 
= sup {r^ + G^{a)} 

?6[0,m) 

= sup {r{ + supjav — ^^(i')}} by (A.7) and Lemma 3.9, 

— sup sup{r^ + av — G^{v)}. (3.61) 

We now argue that sup^gjjlav — G*^(v)} can be achieved when /i^(^) < v < +oo. 

Equation (3.58) gives the values of G*^{v) for v G [h^{^), oo) and we see that G*^{v) 
is differentiable for v e (/i^(^^),oo) values. The derivative (using (3.58) and (3.54)) is 

d„.. J, . *o(As'('f)-0-S'o('!s'(i')) 



The right derivative of G*^ SiS v ^ ^€(^) is 

^^^A , ^o(/^7'(^)-0-^o(/^.'(t^)) 



. ))-\l+ J^^^V,,(;. + e-/.,-^(.))-^0(/. -/.,-(.)) 



^^,^^*o(^)-*o(^ 



*o(^)-*o(^). 
= 0. (3.63) 

Recall that G*^{v) is a convex function of v and note that its subdifferential set (see [24]) 
is a subset of [0, t]. Since a > 0, this implies that the supremum cannot be attained for 

Therefore, 



J(t,t)-ei{r) = sup sup {r^ + av - Gl{v)}. (3.64) 

e6[0,/i) i;e[/i5((M+0/2),oo) 
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To compute Js,t(^) we can assume without loss of generality that s < t (otherwise 
observe that Js,t{'^') = Jt,s{''^))- Then let a — t — s, observe that Js^t^f") = •Jt-{t-s)A^) 
use (3.64). 

□ 

Remcirk 3.10. For s — t — 1, the rate function has a unique at ro — — 2^o(a*/2). 
Even though an explicit Taylor expansion is not easily computable, we can still show the 
first term in the expansion around ro has order 3/2. To simplify the calculations that 
follow, assume that n — 2. 

Let r = ro+e. Using (3.3), it is not hard to verify that the maximizing v is h^{l+^/2). 
Then 

Ji,,{r) = sup {(ro + £)e - 2(logr(l - {/2)) - logr(l + C/2))}. (3.65) 

0<^<2 

Take the ^-derivative of the expression in the braces of (3.65) to show (also by using 
(A.l)^ that the maximizing ^, Cmax, solves the equation 



°^ / 2 2 2 \ 

e = 27-^'o(l- W2)-^o(l+W2) = ^ (-k^ ok-^ + 2A- + f ' 



(3.66) 



All terms in (3.66) are positive, so ^max must be such so that the first term of the series 
satisfies 

^ Smax ^ ' smax 

On the other hand, since ^max < 1; one can easily bound the sum in (3.66) from above, 
to obtain 

^ ^ S ^'^^^A; - l){2k + 1) ^ A {2x - 1)3 " ^"^"^ ^^-^^^ 
Solve (3.67); (3.68) to obtain 

< Uax < 2e'/^ (3.69) 
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[nsj 



[ns\ - [nail 



[nsJ - [nflj+ij 




[nail [nai+i\ [ns\ 



Figure 5: The partition function in the thick-set square bounds from above the product 
of the partition functions of the shaded parts. This is used in the proof of Theorem 3.3. 

We now bound (3.65)m the following manner: 

Jl,l{r) = ro^max + eimax " 2(logr(l - ^,naa/2)) - logr(l + ^„^ax/2)) 

= sup Ui + ei + 2^,{l)i + ^^e + 0{e)]. 

since by (3.69), we can Taylor expand the braces for small values of ^ . Recall that 
To = — 2\l/o(l) and that \l/2 is finite away from 0, to conclude that there exist positive 
constants C\, C2 so that 

Ci(r - rof'^ + 0{{r - r^f'^) < J,^,{r) < C^ir - r^f'^ + 0[[r - r^f'^). (3.70) 
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3.4 Exact free-endpoint rate function 

We now turn our attention to the free endpoint case with no boundary. The conclusion 
of Theorem 3.1 suffices to prove Theorem 3.3. 

Proof of Theorem 3.3. It is easy to check the following bounds for < A; < [^-sj: 

l0g(^Lnij.Ln.J-Lnij)^l°g^Kj^ (3.71) 

< log(ns + 1) + log ( max Zk,\ns\-k) ■ (3.72) 

k 

To get the upper bound let n — > oo using (3.71). 

We show the lower bound. Let £ > and let n large enough so that log(ns + 1) < 

e. 

P{logZ[°*j > nr} < P{logns + log (maxZfc,Ln.j-fe) > nr} 

< P{max logZk,insi-k > n{r - e)} 

k 

lns\ 

< Y,n^ogZk,insi-k > n{r - e)}. (3.73) 

fc=0 

After taking logarithms on both sides of (3.73), we have the upper bound 

logP{logZg*j > nrs} < max logP{logZfe,Ln.j-ik > n{r - e)} + \og{ns + 1). (3.74) 

Now take e' > and a partition = Oi < 02 < ... < = s of the interval [0, s] with 
small enough mesh so that for r fixed, \Jai,s-ai{r) — Jai+i,s-ai{r)\ < s . We conclude that 
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for n large enough (3.72) implies 

n-MogP{logZ[°*j > nr} < 

< max n"MogP{logZL„a,j,Ln8j-Lnaij > n{r - e)} + £ 

[na.i+ij 

< max n"MogP| logZL„a j^Lnsj-Lna^j - yZ log^i,LnsJ-L"aiJ > n{r - s) \ + s. 

i=\nai\+l 

(3.75) 

The same arguments as in (3.39)-(3.43) turn (3.75) into 

lim logP{log > nr} < - rain J„,+„.-„,(r - 2^) + £ 

<- min Ja„s-aAr)+e' (3.76) 

l<i<m 

< - inf J„,,_„(r)+£'. (3.77) 

0<a<s 

Equation (3.76) is a result of letting £ — > and after that adjusting the index of the rate 
function. Finally, let e' — )■ 0, to conclude 

- lim n-MogP{logZ*°* > nr} > inf J„,,_a(r). (3.78) 

n->oo ■' 0<a<s 

By Theorem 2.5, we have convexity of Js,t{^) in the {s,t) argument. This gives 

Js/2,s/2{r) < ya,s-a{r) + | Js-a,a(r) = Ja,s-a{r), (3.79) 

where the last equality follows from the symmetry relation Js,t{'r) = Jt,s{j')- This con- 
cludes the proof. □ 
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Chapter 4 

Multiphase TASEP-Introduction 
and Results 

4.1 Introduction 

The last two chapters study hydrodynamic hmits of totally asymmetric simple exclu- 
sion processes (TASEPs) with spatially inhomogeneous jump rates given by functions 
that are allowed to have discontinuities. We prove a general hydrodynamic limit and 
compute some explicit solutions, even though information about invariant distributions 
is not available. The results come through a variational formula that takes advantage 
of the known behavior of the homogeneous TASEP. This way we are able to get ex- 
plicit formulas, even though the usual scenario in hydrodynamic limits is that explicit 
equations and solutions require explicit computations of expectations under invariant 
distributions. Together with explicit hydrodynamic profiles we can present explicit limit 
shapes for the related last-passage growth models with spatially inhomogeneous rates. 

The class of particle processes we consider are defined by a positive speed function 
c{x) defined for a: e R, lower semicontinuous and assumed to have a discrete set of 
discontinuities. Particles reside at sites of Z, subject to the exclusion rule that admits at 
most one particle at each site. The dynamical rule is that a particle jumps from site i to 
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site i + 1 at rate c{i/n) provided site i + 1 is vacant. Space and time are both scaled by 
the factor n and then we let n — > oo. We prove the almost sure vague convergence of the 
empirical measure to a density p{x,t), assuming that the initial particle configurations 
have a well-defined macroscopic density profile po- 

From known behavior of driven conservative particle systems a natural expectation 
would be that the macroscopic density p{x, t) of this discontinuous TASEP ought to be, 
in some sense, the unique entropy solution of an initial value problem of the type 

Pt + {c{x)f{p)), = 0, p{x, 0) = po{x). (4.1) 

Our proof of the hydrodynamic limit does not lead directly to this scalar conservation 
law. We can make the connection through some recent PDE theory in the special 
case of the two-phase flow where the speed function is piecewise constant with a single 
discontinuity. In this case the discontinuous TASEP chooses the unique entropy solution. 
We would naturally expect TASEP to choose the correct entropy solution in general, 
but we have not investigated the PDE side of things further to justify such a claim. 

The remainder of this introduction reviews briefly some relevant literature and then 
closes with an overview of the contents of these last two chapters. The model and the 
results are presented in Section 4.2. 

Discontinuous scalar conservation laws. The study of scalar conservation laws 

Pt + F{x,p), = Q (4.2) 

whose flux F may admit discontinuities in x has taken off in the last decade. As with 
the multiple weak solutions of even the simplest spatially homogeneous case, a key issue 
is the identification of the unique physically relevant solution by means of a suitable 
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entropy condition. (See Sect. 3.4 of [16] for textbook theory.) Several different entropy 
conditions for the discontinuous case have been proposed, motivated by particular phys- 
ical problems. See for example [1, 3, 4, 10, 15, 18, 21]. Adimurthi, Mishra and Gowda 
[3] discuss how different theories lead to different choices of relevant solution. An inter- 
esting phenomenon is that limits of vanishing higher order effects can lead to distinct 
choices (such as vanishing viscosity vs. vanishing capillarity). 

However, the model we study does not offer more than one choice. In our case the 
graphs of the different fluxes do not intersect as they are all multiples of f{p) — p{l — p). 
In such cases it is expected that all the entropy criteria single out the same solution 
(Remark 4.4 on p. 811 of [2]). By appeal to the theory developed by Adimurthi and 
Gowda [1] we show that the discontinuous TASEP chooses entropy solutions of equation 
(4.1) in the case where c{x) takes two values separated by a single discontinuity 

Our approach to the hydrodynamic limit goes via the interface process whose limit 
is a Hamilton- Jacobi equation. Hamilton-Jacobi equations with discontinuous spatial 
dependence have been studied by Ostrov [21] via mollification. 

Hydrodynamic limits for spatially inhomogeneous, driven conservative peir- 
ticle systems. Hydrodynamic limits for the case where the speed function possesses 
some degree of smoothness were proved over a decade ago by Covert and Rezakhan- 
lou [14] and Bahadoran [5]. For the case where the speed function is continuous, a 
hydrodynamic limit was proven by Rezakhanlou in [23] by the method of [27]. 

The most relevant and interesting predecessor to our work is the study of Chen et 
al. [10]. They combine an existence proof of entropy solutions for (4.2) under certain 
technical hypotheses on F with a hydrodynamic limit for an attractive zero-range process 
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(ZRP) with discontinuous speed function. The hydrodynamic hmit is proved through a 
compactness argument for approximate solutions that utihzes measure- valued solutions. 
The approach follows [5, 14] by establishing a microscopic entropy inequality which 
under the limit turns into a macroscopic entropy inequality. 

The scope of [10] and our work are significantly different. Our flux F{x,p) = 
c{x)p{l — p) does not satisfy the hypotheses of [10]. Even with spatial inhomogeneities, 
a ZRP has product-form invariant distributions that can be readily written down and 
computed with. This is a key distinction in comparison with exclusion processes. The 
microscopic entropy inequality in [10] is derived by a coupling with a stationary process. 

Finally, let us emphasize the distinction between the present work and some other 
hydrodynamic limits that feature spatial inhomogeneities. Random rates (as for example 
in [27]) lead to homogenization (averaging) and the macroscopic flux does not depend 
on the spatial variable. Somewhat similar but still fundamentally different is TASEP 
with a slow bond. In this model jumps across bond (0, 1) occur at rate c < 1 while 
all other jump rates are 1. The deep question is whether the slow bond disturbs the 
hydrodynamic profile for all c < 1. V. Beffara, V. Sidoravicius and M. E. Vares have 
announced a resolution of this question in the affirmative. Then the hydrodynamic limit 
can be derived in the same way as in the main theorem here. The solution is not entirely 
explicit, however: one unknown constant remains that quantifies the effect of the slow 
bond (see [28]). [6] generalizes the hydrodynamic limit of [28] to a broad class of driven 
particle systems with a microscopic blockage. 

Organization. Section 4.2 contains the main results for the inhomogeneous corner 
growth model and TASEP. Sections 5.1 and 5.2 prove the limits. Section 5.3 outlines the 
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explicit computation of density profiles for the two-phase TASEP. Section 5.4 discusses 
the connection with PDE theory. 

Notational conventions. The Exp(c) distribution has density f{x) = ce~^ for < 
X < oo. Two last passage time models appear: the corner growth model whose last- 
passage times are denoted by G, and the equivalent wedge growth model with last- 
passage times T. H{x) — l[o,oo)(2^) is the Heavyside function. C is a constant that may 
change from line to line. 



The corner growth model connected with TASEP has been a central object of study in 
this area since the seminal 1981 paper of Rost [25]. So let us begin with an explicit 
description of the limit shape for a two-phase corner growth model with a boundary 
along the diagonal. Put independent exponential random variables {i^j^eN^ on the 
points of the lattice with distributions 



4.2 Results 



Exp(ci), 




Exp(ciAc2), ifi=j 



(4.3) 



^ Exp(c2), 



if i > j. 



We assume that the rates satisfy ci > C2. 



Define the last passage time 




max 

7r6n(m,n) 




(4.4) 
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Figure 6: Optimal macroscopic paths that give the last passage time constant described 
in Theorem 4.1. 



where Il{m,n) is the collection of weakly increasing nearest-neighbor paths in the rect- 
angle [m] X [n] that start from (1, 1) and go up to {m,n). That is, elements of n(m,n) 
are sequences {(1, 1) = f i, V2, ■ ■ ■ , Vm+n-i = ijn, n)} such that Vi+i — Vi = (1, 0) or (0, 1). 

Theorem 4.1. Let the rates ci > C2 > 0. Define c = Ci/c2 > 1 and 6 = 2c — 1 — 
2a/c(c — 1). Then the a.s. limit 



^{x,y) = lim n ^G{\nx\, \ ny\) 



exists for all {x,y) G (0,oo)^ and is given by 



ifO<x< iP'y 



^ix.,y^=\, ^--i^ + ^f I {l + hf-^ch- ^fb-y<x<y 

C2^(v/^+v^)\ if y <x <+oo. 

This theorem will be obtained as a side result of the development in Section 5.1. 
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We turn to the general hydrodynamic limit. The variational description needs the 
following ingredients. Define the wedge 

W = {ix,y)eR^:y>0,x>-y} 

and on W the last-passage function of homogeneous TASEP by 

7(a;, y) = {y/x + y + ^f. (4.5) 

Let x(s) = {xi{s),X2{s)) denote a path in and set 

H{x, = {x e C([0, 1], W) : X is piecewise C\ x(0) = (0, 0), 

x(l) = {x,y), x'(s) e yV wherever the derivative is defined}. 

The speed function c of our system is by assumption a positive lower semicontinuous 
function on M. We assume that at each x e M 

c{x) — min \ lim c(y), lim c{y) >. (4.6) 

In particular we assume that the limits in (4.6) exist. We also assume that c{x) has only 
finitely many discontinuities in any compact set, hence it is bounded away from in any 
compact set. 

For the hydrodynamic hmit consider a sequence of exclusion processes r^" = {Vi{t) '■ 
i E t & M_i_) indexed by n G N. These processes arc constructed on a common 
probability space that supports the initial configurations {r^"(0)} and the Poisson clocks 
of each process. As always, the clocks of process 77" are independent of its initial state 
77"^ (0). The joint distributions across the index n are immaterial, except for the assumed 
initial law of large numbers (4.12) below. In the process 77" a particle at site i attempts 
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a jump to i + 1 with rate c{i/n). Thus the generator of 77" is 

Lnfiv) = Yl c{xn-'Hx)(l - v(x + l))(/(77^'^+^) - fiv)) (4-7) 

for cyhnder functions / on the state space {0, 1}^. The usual notation is that particle 
configurations are denoted by 77 = {ri{i) : i e Z) e {0, 1}^ and 

when i = x 
= <J 1 wheni = x + l 
r\{i) wheni7^x,x + l 

is the configuration that results from moving a particle from x to x + 1. Let denote 
the number of particles that have made the jump from site i to site i + 1 in time interval 
[0,t] in the process rp . 

An initial macroscopic profile po is a measurable function on R such that < p^ix) < 
1 for all real x, with antiderivative Vq satisfying 



^0(0) = 0, vo{b)-vo{a)= / po{x)dx. 

J a 

The macroscopic flux function of the constant rate 1 TASEP is 



(4i 



Its Legendre conjugate 



p(l-p), ifO<p<l 
/(P) = <! (4.9) 
—00, otherwise. 



r(y)= inf {yp-/(p)} 

0<p<l 



represents the limit shape in the wedge. We orient our model so that growth in the 
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wedge proceeds upward, and so we use g{y) — —f*{y)- It is explicitly given by 

-y, iiy<-l 
g{y) = sup {/(p) -yp}^{ \{l - y)\ if - l<y<l (4.10) 

0<p<l ^ 

^0, ify>l. 

For a; e M define v{x, 0) = ^0(2;), and for t > 0, 

v{x,t)=SM^[vo{w{Q))- I c(w;(s))(7(^^|^^j rfsj (4.11) 



w{-) I, JO 



where the supremum is taken over continuous piecewise paths w : [0, t] — > R that 
satisfy w{t) = x. The function v{x, t) is Lipschitz continuous jointly in (x, t). (see Section 
5.2) and it has a derivative almost everywhere. The macroscopic density is defined by 

p{x,t) = v^{x,t). 

The initial distributions of the processes 77" are arbitrary subject to the condition 
that the following strong law of large numbers holds at time t — 0: for all real a <b 

lim - y 77^(0) = / po(x)dx a.s. (4.12) 

i=[na}+l 

The second theorem gives the hydrodynamic limit of current and particle density for 
TASEP with discontinuous jump rates. 

Theorem 4.2. Let c{x) he a lower semicontinuous positive function satisfying (4.6), 

with finitely many discontinuities in any compact set. Under assumption (4.12), these 
strong laws of large numbers hold at each t > 0: for all real numbers a < b 

lim Jll^^Ant) = Vo{a) —v{a,t) a.s. (4-13) 

and 

lim- rflint)^ / p{x,t)dx a.s. (4-14) 

i=[na}+l 
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where v{x,t) is defined by (4.11) and p{x,t) — Vx{x,t). 

Remcirk 4.3. In a totally asymmetric K -exclusion with speed function c the state space 
would be {0, 1, ... , K}^ with K particles allowed at each site, and one particle moved 
from site x to x + 1 at rate c{x/n) whenever such a move can be legitimately completed. 
Theorem 4-2 can be proved for this process with the same method of proof . The definition 
of the limit (4.11) would be the same, except that the explicit flux f and wedge shape g 
would be replaced by the unknown functions f and g whose existence was proved in [27]. 

To illustrate Theorem 4.2 we compute the macroscopic density profiles p{x, t) from 
constant initial conditions in the two-phase model with speed function 

c{x)=ci{l-H{x)) + C2H{x) (4.15) 

where H{x) = l[o,oo)(a;) is the Heavyside function and Ci > C2. (The case Ci < C2 can 
then be deduced from particle-hole duality.) The particles hit the region of lower speed 
as they pass the origin from left to right. Depending on the initial density p, we see the 
system adjust to this discontinuity in different ways to match the actual throughput of 
particles on either side of the origin. The maximal flux on the right is C2/4 which is 
realized on the left at densities p* and I — p* with 

P* = I - |\/l -C2/C1. 

Corollary 4.4. Let Ci > C2 and the speed function as in (4.15). Then the macroscopic 
density profiles with initial conditions po{x, 0) — p are given as follows. 
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1 

2 










P 











C2{l-2r*)t 



C2{l-2p)t 



Figure 7: Density profile p{x,t) in the two-phase (ci > C2) TASEP when we start from 
constant initial configurations po{x) = p € (0,p*). 



(i) Suppose < p < p*. Define r* = r*{p) = \- - 4p(l - p)ci/c2- Then 



if — 00 < X < 

z/0 < X < C2(l - 2r*)t 

if €2(1 - 2r*)t <x < C2(l - 2p)t 

if (1 — 2p)tc2 < X < +00 



p(x,t) 



1 



X 



2 V ^ tC2 



p 



(ii) Suppose p* < p < \. Then 

P 



p{x,t) = < 



1-p* 
1 



X 



tC2 



if — 00 < X < — tci(p — p* 
if -tei(p-p*) < X < 
z/0 < X < (1 - 2p)tC2 

if (1 - 2p)tc2 < X < +00 



(iii) Suppose p > h Define r* = r*{p) 



1 _ 1 

2 2 



1 — 4p(l — p)c2/ci. Then 



p if — 00 < X < — tci(p — r*) 

p(x,t)=<j l-r* if -tci{p-r*) <x <0 
p if < X < +00 



(4.16) 



(4.17) 



(4.18) 
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1-p* 
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. 2_ _ 






P ^^^^^ 









ci{p* ~ p)t C2(l-2p)t 

Figure 8: Density profile p{x,t) in the two-phase (ci > C2) TASEP when we start from 
constant initial configurations po{x) = p G [p*, |]. 





1-r* 




p 








2 





ciij-* - p)t 

Figure 9: Density profile p{x,t) in the two-phase (ci > 02) TASEP when we start from 
constant initial configurations po{x) = p G !)• 

Remark 4.5. Taking t 00 in the three cases of Corollary 4-4 gives a family of 
macroscopic profiles that are fixed by the time evolution. A natural question to investigate 
would be the existence and uniqueness of invariant distributions that correspond to these 
macroscopic profiles. 

Next we relate the density profiles picked by the discontinuous TASEP to entropy 
conditions for scalar conservation laws with discontinuous fluxes. The entropy conditions 
defined by Adimurthi and Gowda [1] are particularly suited to our needs. Their results 
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give uniqueness of the solution for the scalar conservation law 

Pt+{F{x,p))^^0, xeR,t>0 

p{x, 0) = po{x), x e M 
with distinct fluxes on the half-lines: 

F{x,p) = H{x)Mp) + (1 - H{x))Mp) 



(4.19) 



(4.20) 



where fr,fi G C^(]R) are strictly concave with superlinear decay to — oo as — >■ oo. 
A solution of (4.19) means a weak solution, that is, p e I/^c(IR x M+) such that for all 
continuously differentiable, compactly supported test functions e (R x R+), 

f°^£ + ^^^^^^ + j^"" p{x,0)cf>(x,0)dx^O. (4.21) 

(4.21) is the weak formulation of the problem 

Pt + fr{p)x = 0, iorx>0,t>0 



(4.22) 



Pt + fe{p)x = 0, forx<0,i>0 
fr{p{0+,t)) = fe{p{0-,t)) for a.e. t > 
p{x,0) = po{x). 

The entropy conditions used in [1] come in two sets and assume the existence of 
certain one-sided limits: 

(Ei) Interior entropy condition, or Lax-Oleinik entropy condition: 



p{x+, t) > p{x—, t) for X 7^ and for all i > 0. 



(4.23) 
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(Ef,) Boundary entropy condition at x — 0: for almost every t, the hmits p{0±,t) 
exist and one of the following holds: 



/;(p(o+,o)>o 


and 


/;(p(o- 


>o, 


(4.24) 


/;(p(o+,t))<o 


and 


/;(p(o- 


<o, 


(4.25) 


/;(p(o+,i))<o 


and 


/;(p(o- 


>0. 


(4.26) 



Define 

G,{p) = i{x > o}/;(p) + i{x < o}/;(p) + i{x = o}min(/;(o),/;(o)), 

where /* and // are the convex duals of fr and fi. Set Vq{x) = Jq Po{0) dO and define 

V{x,t)=sv.^{vo{w{Q))+ f G^^s){w\s))ds\ (4.27) 
w{-) K JO ) 

where the supremum is over continuous, piecewise linear paths w : [0, — > R with 

w{t) = X. 

Theorem 4.6. [1] Let po G L°^(]R) and define V by (4.27). Then V is a uniformly 
Lipschitz continuous function and p{x, t) — Vx{x, t) is the unique weak solution of (4.22) 
that satisfies the entropy assumptions (Ei) and {Eb) in the class L°° fl BV\oc and with 
discontinuities given by a discrete set of Lipschitz curves. 

It is easy to check that the two-phase density profile p{x, t) in Corollary 4.4 is a weak 
solution (in the sense of (4.21)) to the scalar conservation law (4.19) with flux function 
F{x,p) = c{x)p{l — p). However we cannot immediately apply this theorem in our case 
since the two-phase flux function F(a;, p) = (1 — H[x))cif{p) + H{x)c2f{p) is flnite only 
for p e [0,1] and in particular is not C^. We show how we can replace F{x,p) with 
F{x,p) in the above theorems in Section 5.4. In particular, we prove the following. 
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Theorem 4.7. For p e R define fr{p) — 02(1 — p)p and /^(p) = Ci(l — p)p to be the flux 
functions for the scalar conservation law (4.22). Let the initial macroscopic profile for 
the hydrodynamic limit be a measurable function < po{x) <1. Then the macroscopic 
density profile p{x, t) from the hydrodynamic limit in Theorem 4-2 is the unique solution 
described in Theorem 4-6. 
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Chapter 5 

Multi-phase TASEP 

5.1 Wedge last passage time 

The strategy of the proof of the hydrodynamic hmit is the one from [27] and [28] . Instead 
of the particle process we work with the height process. The hmit is first proved for 
the jam initial condition of TASEP (also called step initial condition) which for the 
height process is an initial wedge shape. This process can be equivalently represented 
by the wedge last-passage model. Subadditivity gives the limit. The general case then 
follows from an envelope property that also leads to the variational representation of 
the limiting height profile. In this section we treat the wedge case, and the next section 
puts it all together. 

Recall the notation and conventions introduced in the previous section. In partic- 
ular, c{x) is a positive, lower semicontinuous speed function with only finitely many 
discontinuities in any compact set. Define a lattice analogue of the wedge W by 

C = {{t,j)eZ':j>l,t>-j + l} (5.1) 

with boundary dC = {(i, 0) : i > 0} U {(i, -i) : i < 0}. 

For each n e N construct a last-passage growth model on C that represents the 
TASEP height function in the wedge. Let {Tj"- : e >C}ngN denote a sequence of 

independent collections of i.i.d. exponential rate 1 random variables. We need an extra 
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index I to denote the shifting. Define weights 

= \ (^i)e^- (5-2) 

For £ e Z and n e N assign to site e C the random variable uj^frPj. Given 

lattice points (a, 6), {u,v) G £, n((a,6), {u,v)) is the set of lattice paths n = {{a,b) = 
(io, Jo)) {ip-ijp) — {'^■i'^)} whose admissible steps satisfy 

{iuji) - (ii-iJi-i) e {(1,0), (0, 1), (-1, 1)}. (5.3) 

In the case that (a, 6) = (0,1) we simply denote this set by Il{u,v). For {u,v) e £, 
£ e M and n G N denote the wedge last passage time 

T-'\u,v)= max uj^r^:, (5.4) 

with boundary conditions 

T"'^(k,i;) = for {u,v)edC. (5.5) 

Admissible steps (5.3) come from the properties of the TASEP height function. No- 
tice that (0, 1) is in fact never used in a maximizing path. 

To describe macroscopic last passage times define, for (x, y) G W and g G M, 

n(a:,y)= sup [ f -^l^^p^ ds]. (5.6) 
.r(.\czv(^ „\ ^ In C{X^{S] — a] J 



^{■)&H{x,y) ^ Jo c{Xi{s) - q) 

Theorem 5.1. For all q and {x,y) in the interior ofW 



lim n-^T'^'L'*^J([nxJ, [nyj) = r'i{x,y) a.s. (5.7) 
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Remcirk 5.2. In a constant rate c environment the wedge last passage limit is 

lim V(M,LnyJ) = c-S(x,y) = c-i(v/^+v/^)'. (5.8) 

n— )-oo fi 

The limit 7(x,y) is concave, hut this is not true in general for V^{x,y). In some special 
cases concavity still holds, such as if the function c{x) is nonincreasing if x < and 
nondecreasing if x > 0. 

To prove Theorem 5.1 we approximate c{x) with step functions. Let — oo = Oi < 
a2 < ■■■ < ttL^i < — +00, and consider the lower semicontinuous step function 

L-1 L-1 
m=l m=2 

Proposition 5.3. Let c{x) be given by (5.9). Then limit (5.7) holds. 

On the way to Proposition 5.3 we state prehminary lemmas that will be used for 
pieces of paths. We write Cj for the rate values instead of to be consistent with the 
notation in Theorem 4.1. 

Lemma 5.4. Assume that there is a unique discontinuity a2 — for the speed function 
c{x) in (5.9). Then for y > 

lim n-'T-'\0, lny\) = . . = ds a.s. 

n-^oo min{ci,C2j Jo c(0) 

Proof. The upper bound in the limit is immediate from domination with constant rates 
c(0). 

For the lower bound we spell out the details for the case ci > C2. Let £ > 0. To 
bound T"'°(0, lny\) from below force the path to go through points (0, 1), {([^t-I/^J, {k — 
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[eyn\ 2 [eyn\ 



Figure 10: A possible microscopic path forced to go through opposite corners of the 
parallelograms R^. The striped area left of x = is the Ci-rate region. 



1) [nye\ ) : A; = 1, . . . , [e'^J } and (0, [ny\). For 1 < A; < [e'^] let r"(i?^) be the last pas- 
sage time from ( [nye\ ,{k — 1) [nye\ ) to ( [nye\ , k [nye\). RJ^ refers to the parallelogram 
that contains all the admissible paths from {Ynye\,{k — l)Ynye\) to {\nye\,k\nye\). 
Each lies to the right of x = and therefore in the C2-rate area. (See Fig. 10.) 

Let < 5 < £:C2~^7(0, y). A large deviation estimate (Theorem 4.1 in [26]) gives a 
constant C = C{c2,y,e,6) such that 



F{T:^iRl)<niec,-'^iO,y)-6)}<e 



(5.10) 



By a Borel-Cantelli argument, for large n, 

Le-ij-i 

T^'°(0, [ny\) > J2 T^'iRl) > n{[e~'\ - l){ec,~HQ,y) - 6). 

k=l 

This suffices for the conclusion. 



□ 
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Remcirk 5.5. This lemma shows why it is convenient to use a lower semi- continuous 
speed function. A path that starts and ends at the same discontinuity stays mostly in the 
low rate region to maximize its weight. This translates macroscopically to the formula 
for the limiting time constant obtained in the lemma, involving only the value of c at the 
discontinuity. If the speed function is not lower semi- continuous, we can state the same 
result using left and right limits. 

Lemma 5.6. Let a — < b < +oo be discontinuities for the step speed function c{x) 
and c{x) — r for a < x < b. Take z e [0, b] . Let T"'{\nz\ , \ ny\ ) be the wedge last passage 
time from (0, 1) to {lnz\, lny\) subject to the constraint that the path has to stay in the 
r-rate region (a, b) x (0, +oo), except possibly for the initial and final steps. Then 

lim n-^f"([nzj, [ny\) = r"^7(z,y) a.s. (5.11) 

n— )-oo 

Same statement holds if b < z < 0. 

Proof. The upper bound limn~^T'"'([n2;J, [ny\) < r~^j{z,y) is immediate by putting 
constant rates r everywhere and dropping the restrictions on the path. For the lower 
bound adapt the steps of the proof in Lemma 5.4. □ 

Lemma 5.6 is a place where we cannot allow accumulation of discontinuities for the 
speed function. 

Before proceeding to the proof of Proposition 5.3 we make a simple but important 
observation about the macroscopic paths x(s) = {x^{s),x^{s)), s e [0, 1], in 1-L{x,y) for 
the case where c{x) is a step function (5.9). 

Lemma 5.7. There exists a constant C — C{x,y,c{- — q)) such that the supremum in 
(5.6) comes from paths in 'H(a;,7/) that consist of at most C line segments. Apart from 
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Figure 11: A possible microscopic path described in Lemma 5.6. The path has to stay 
in the unshaded region. 

the first and last segment, these segments can be of two types: segments that go from 
one discontinuity of c{- — q) to a neighboring discontinuity, and vertical segments along 
a discontinuity. 

Proof. Path x is a union of subpaths {xj} along which c(xj(s) — q) is constant, except 
possibly at the endpoints. Given such a subpath (xj(s) : tj < s < tj+i), concavity of 7 
and Jensen's inequality imply that the line segment (pj that connects ^j(tj) to Xj(tj+i) 
dominates: 



Consequently we can restrict to paths that are unions of line segments. 

To bound the number of line segments, observe first that the number of segments 
that go from one discontinuity to a neighboring discontinuity is bounded. The reason 
is that the restriction x'(s) G W forces such a segment to increase at least one of the 
coordinates by the distance between the discontinuities. 

Additionally there can be subpaths that touch the same discontinuity more than once 
without touching a different discontinuity. Lower semi-continuity of c(-) and Jensen's 
inequality show again that the vertical line segment that stays on the discontinuity 
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dominates such a subpath. Consequently there can be at most one (vertical) line segment 
between two line segments that connect distinct discontinuities. □ 

Next a lemma about the continuity of F*. We write r^((a, 6), {x,y)) for the value in 
(5.6) when the paths go from (a, b) to {x, y) e (a, h) + W. 

Lemma 5.8. Fix z,w > Q. Then there exists a constant C = C{z, w, c(- — q))<oo such 
that for all < S < 1 and < a < z 

n((a, 0), {z, 5)) - n((a, 0), {z, 0)) < CV~5, (5.12) 

and for < b < w 

n((-6, 6), {-w, w + S))- r«((-6, 6), {-w, w)) < CV5. (5.13) 



Proof. Pick 6 e (0, 1] and consider the point {z,S) in W. For any x = {x^{s),x'^{s)) e 
^{z, 6) set 

7(x,g)= /Iff , ds. (5.14) 
Jo c{x^{s)-q) 

Let £ > and assume that 4> — {(p^,4>'^) G ^(^^ <^) is a path such that r^(^, 5) — 1{4>, q) < 
e. Lemma 5.7 implies that we can decompose into disjoint linear segments 0j so that 
4> = "^jLi 4>j and 4>j : [sj-i, Sj] W. Here (pj means path concatenation. 
We can find segments (f)j(k): ^ ^ k < N, such that 

(l^](k)(Sj{k)-l) < 0](ifc)(Sj(fc)), 0](fe)(Sj(fe)) = (t>](k+l)(Sj{k+l)-l), 

0j^^^(sj(i)_i) = 0, and 0](7v)(*j(iv)) — ^- In other words, the projections of the segments 
4>j{k) cover the interval [0, z\ without overlap and without backtracking. 
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We bound the contribution of the remaining path segments to I{(l),q). Let J = 
U'l^'-^lsj^k), Sj(k+i)-i] be the leftover portion of the time interval [0,1]. The subpath 
4>{s), s e [sj(fe), Sj(fe+i)_i], (possibly) eliminated from between and 4>j{k+i) satisfies 
(f)^{sj(^k)) = Note that 7(0, 6) < 2a+46for {a,b) e Wand /o^(0^)'(s) ds = 5. 

We can bound as follows: 

<C [ 2{(j)^y{s)ds + C [ 4{(j)^y{s)ds 
J J Jo 

= + 4C5. (5.15) 

Set tk — Sj(fe)-i < life = SjXfe)- Define a horizontal path w from (0,0) to (-2,0) with 
segments 

Wk{s) = (0](fe)(s),O), for tk<s< Uk, (5.16) 
and constant on the complementary time set J. 
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To get the lemma, we estimate 

J J c{<t>\s) - q) J[o,i]\j c{(f)\s) - q) 

N 

<cs + J2 {H<Pj{khQ) - n^k, q)) + r^iz, 0) 

k=l 

<C5 + C'J2 Wm(s)) - i(w',{s))) ds + 0) 
k=i -^^k 

k=i 

^ ^Jmmi^)^i<l>%k)i■^) + {<l>%,^{s)ds^+T'^{z,0) 
<CS + C'J2 [[\<P%(k)is) ds^ m'm^^) + (^)) ds^ ' + T\z, 0) 



+ 2 

AT 



<C5 + C'(J2 i^%k)is)ds) X 
k=i -^^k 



V 1.-1 Jtu / 



k=i 



fe=i '^*'= 

{{<P%ik){s) + {<P')',(U){s))ds + T\z,Q) 

<C5 + C'V8\/T+5 + T\z, 0) 
< C5 + C"v^v^ + n(^,0). 

The first inequahty (5.12) follows for a = by letting £ go to 0. It also follows for all 
a e [0, by shifting the origin to a which replaces z with z — a. 

For the second inequality (5.13) the arguments are analogous, so we omit them. □ 

Corollary 5.9. Fix {x,y) e W. Then there exists C = C{x,y,c{- — q)) < oo such that 
for allO<S<l 

Ti{x, y + 5)- n(a;, y) < CVS. (5.17) 
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(0,0) (z,0) {z + 6,0) 

Figure 12: A possible macroscopic path from (0,0) to {z,6). The dotted vertical lines 
are discontinuity columns of c(- — q). The error from eliminating the segments outside 
the two vertical dashed lines and from eliminating pathologies (like the circled part) is 
of order 6 and a comparison with the horizontal path leads to an error of order \/6. 



Proof. Let A{{a,b), {x,y)) be the parallelogram with sides parallel to the boundaries 
of the wedge, north-east corner the point {x,y) and south-west corner at {a,b). If 
(a, 6) = (0,0) we simply write A{x,y). 

Let e > 0. Let 0^ a path such that V^Xjy + 6) — 1(0'^, g) < e. Let u be the point 
where 0^ first intersects the north or the east boundary of A{x,y). Without loss of 
generality assume it is the north boundary and so u = (a, y) for some a G [—y, x]. Then, 

T%x,y + S)-e<I{f,q) 

<n(a,y) + n((a,2/),(x,y + 5)) 

= n(a, y) + n((a, y), {x, y)) + n((a, y), (x, y + 5)) - r^((a, y), (x, y)) 
< V\x, y) + V\{a, y), (x, y + 6)) - ^((a, y), (x, y)). (5.18) 

The last inequality gives 

n(x, y + 6)- n(x, y) < n((a, y), (x, y + 6)) ~ n((a, y), (x, y))+e<CV6 + e (5.19) 

by Lemma 5.8. Let e decrease to to prove the Corollary. □ 



93 

Proof of Proposition 5.3. Fix {x,y) in the interior of W. For x = {x^{s),x^{s)) e 
T-L{x,y) set 

7(x,g)= Tv^^ci.. (5.20) 



c{x^{s) - q) 



We prove first 



lim n-^T'^'L'^^J(Lna;J, LnyJ) > n(a;,y) = sup /(x,g). (5.21) 
It suffices to consider macroscopic paths of the type 

H 

x(s) = E^^(^)lfe.^.+i)(«) (5-22) 

where H & N, x^ is the straight hne segment from x(sj) to x(sj+i), c{xi{s) — q) — rm^ 
is constant for s G (sj,Sj+i), and by continuity Xj(sj_|_i) = Xj+i(sj+i). 

Let tt" be the microscopic path through points (0, 1), {\n'^j{sj)\ : 1 < j < H} and 
( [nx\ , [ny\ ) constructed so that its segments tt" satisfy these requirements: 

(i) TTj lies inside the region where (jJ^'}"''^^ = r^. is constant, except possibly for the 
initial and final step; 

(ii) tt" maximizes passage time between its endpoints \n'^j{sj)\ and Lnxj+i(sj+i)J 
subject to the above requirement. 

Let 

77'L"^J=m^ax (5-23) 

denote the last-passage time of a segment subject to these constraints. Observe that the 
proofs of Lemmas 5.4 and 5.6 do not depend on the shift parameter q, therefore 

hm n-^T^'L^^J = 7(x,(g,)-x,+i(g,)) ^ li^'jjs)) 



r-mj Jsj c{x]{s)-q) 
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Adding up the segments gives the lower bound: 

lim n-^r"'L"«J([nxJ, [ny\) > lim J^^'^T^'^^'^ 



n— )-oo 

J 



•/ s. 



^^■+1 7(x;(s)) 

ds — I{x,q). 



c{x^{s) — q) 
Now for the complementary upper bound 

]hR n-^r"'L"«J([nxJ, lny\) < T'^(x,y). (5.24) 

n— >oo 

Each microscopic path to ([nxj, [ny\) is contained in nA for a fixed parallelogram 
A C. yv with sides parallel to the wedge boundaries. Pick e > 0. Let > be a 
lower bound on all the rate values that appear in the set A. Find S > such that 
\j{v) — ^{w)\ < £r* for all v,w E A with \v — w\ < S and 5 < 1 so that Corollary 5.9 is 
valid. 

Consider an arbitrary microscopic path from (0, 1) to ([nxj, [riyj). Given the speed 
function and q, there is a fixed upper bound Q = Q{x,y) on the number of segments 
of the path that start at one discontinuity column ([nojj + [ngj) x N and end at a 
neighboring discontinuity column ( [nojiij + \nq\ ) x N. The reason is that there is an 
order n lower bound on the number of lattice steps it takes to travel between distinct 
discontinuities in nA. 

Fix K eN and partition the interval [0, y] evenly by bj = jy/K, < j < K, so that 
y/K < S/Q. Make the partition finer by adding the y— coordinates of the intersection 
points of discontinuity lines {a, + q} x 1R+ with the boundary of A. 

Let tt" be the maximizing microscopic path. We decompose tt" into path segments 
{tt" : < j < Mn} by looking at visits to discontinuity columns ([najj + [nq\) x N, both 
repeated visits to the same column and visits to a column different from the previous 
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one. Let {0 = b^^ < h^^ < hk^ < ... < h^^n-i — ^kM„ = 2/} be a sequence of partition 
points and {0 = Xq, Xi = a^^ + q, X2 = + g, . . . ,xm„ = x} a sequence where Xj 
for < J < Mn are discontinuity points of the shifted speed function c(- — q). We can 
create the path segments and these sequences with the property that segment tt" starts 
at {lnxj\,l) with / in the range [nfe^J < / < [nhk^+il and ends at {lnxj+i\,l') with 
L^^ifcj+i J < ^' < b^^kj+i+il ■ In an extreme case the entire path tt" can be a single segment 
that does not touch discontinuity columns. 

In order to have a fixed upper bound on the total number of segments, uniformly 
in n, we insist that for < j < Mn — 1 the labels satisfy: 

(i) For odd j, tTj starts and ends at the same discontinuity column ([nx^J, • ). The 
rate relevant for segment tt" is r^. — c{amj)- 

(ii) For even j, tt" starts and ends at different neighboring discontinuity columns, 
and except for the initial and final points, does not touch any discontinuity column and 
visits only points that are in a region of constant rate ri. . 

The above conditions may create empty segments. This is not harmful. Replace Q 
with 2Q + 2 to continue having the uniform upper bound M„ < Q. 

Let T{'k'^) be the total weight of segment tt^ . Let tt" be the maximal path from 
([nxjj, \ nbkj\) to ([nxj+ij, L^^'fe^+i+iJ) in an environment with constant weights uiij — 
r^^ everywhere on the lattice, with total weight TJ'. T-' > T{-kJ), up to an error from 
the endpoints of tt^. 

Theorem 4.2 in [26] gives a large deviation bound for TJ*. Consider a constant rate 

r environment and the maximal weight T(([nMiJ, [nfij), (LnM2j, [^"^2])) between two 
points {ui,Vi) and {u2, V2) such that their lattice versions can be connected by admissible 
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paths for all n. Then there exists a positive constant C such that for n large enough, 

P|r(([niiij, [nvi\),{[nu2\,[nv2\)) > nr-'^-f{u2 - Ui,V2 - Vi) + ne^ < e"^". (5.25) 

There is a fixed finite collection out of which we pick the pairs {{xj, bk^), (3:^^+1, ^fe^+i+i)} 
that determine the segments tt^. By (5.25) and the Borel-Cantelli lemma, a.s. for large 
enough n, 

77 < nrj^-f{xj+i - Xj, bkj^^+i - 6fc.) + ne for < j < M„. (5.26) 

Define 5i > by y + 5i = Y.f=o~\hj+i+i - hj)- Since y = E^'o'^l^fc.+i " ^fcj 
by the choice of the mesh of the partition {6^}, we have 5i < Mn5/Q < 5. Think of 
{xjjf-i — Xj, — bk-), < j < M„, as the successive segments of a macroscopic path 

from (0, 0) to (x, y + Si). 

For sufficiently large n so that (5.26) is in effect, 

r"'L-''J ( [nx\ , [ny\ ) < J] ^7 < n r^^lixj+i - Xj, - + nQe 

i=i i=i 

< nr'i{x,y + Si)+nQe 

< nn(a;, y) + nCVS + nQe. 

The last inequality came from Corollary 5.9. Let S ^ 0. Since e was arbitrary the upper 
bound (5.24) holds. □ 

Proof of Theorem 5.1. Fix {x,y). For each £ > we can find lower semicontinuous step 
functions ci and C2 such that ||ci — C2II00 < £ and on some compact interval, large enough 
to contain all the rates that can potentially influence r'^{x,y), ci{x) < c{x) < C2{x). 
When the weights in (5.2) come from speed function q let us write Tj for last passage 
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times and Fj for their limits. An obvious coupling using common exponential variables 
gives 

Tr'L"^J(M, M) > M) > lny\). 

Letting a > denote a lower bound for c{x) in the compact interval relevant for {x,y), 
we have this bound for x G l-i{x,y): 

Jo ^ci[xi[s)-q) C2(xi(s) - g) J Jo cf(xi(s)-g) 

< £Q;~^7(x,y). 

Therefore the limits also have the bound 

Q <V\{x,y) -Vl{x,y) <C{x,y)e. 

From these approximations and the limits for % in Proposition 5.3 we can deduce The- 
orem 5.1. □ 

Proof of Theorem 4-1 ■ We can construct the last passage times G{x,y) of the corner 
growth model (4.4) with the same ingredients as the wedge last passage times T"'°(a;, y) 
of (5.4), by taking Y(jj) = ^7-j,j'^i'-3,3- Then T'*'°(x, y) = G{x + y, y) and we can transfer 
the problem to the wedge. The correct speed function to use is now c{x) — Cil{x < 
0} + C2l{x > 0}. In this case the limit in Theorem 5.1 can be solved explicitly with 
calculus. We omit the details. □ 



5.2 Hydrodynamic limit 

In this section we sketch the proof of the main result Theorem 4.2. This argument is 
from [27, 28]. 
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5.2.1 Construction of the process and the variational coupHng 

For each n e N we construct a Z- valued height process z'^{t) — {zf{t) : i e Z). The 
height values obey the constraint 

< ^r+iW - z^t) < 1. (5.27) 

Let {Vf} be a collection of mutually independent (in i and n) Poisson processes with 
rates c" given by 

= c{n~H), (5.28) 

where c{x) is the lower semicontinuous speed function. Dynamically, for each n and i, 
the height value is decreased by 1 at event times of P", provided the new configuration 
does not violate (5.27). 

After we construct z'"{t), we can define the exclusion process 77" (t) by 

V?{t)^zm-^tiit)- (5-29) 
A decrease in z^ is associated with an exclusion particle jump from site i to i + 1. Thus 
the process keeps track of the current of the r^^-process, precisely speaking 

Assume that the processes 2;" have been constructed on a probability space that 
supports the initial configurations z'^{0) — (-^"(0)) and the Poisson processes {Df} that 
are independent of (zf (0)). Next we state the envelope property that is the key tool 
for the proof of the hydrodynamic limit. Define a family of auxiliary height processes 
l^n.fc : n eN, k e Z} that grow upward from wedge-shaped initial conditions 

0, if i > 

C'(0) = <! (5.31) 
-i, if i < 0. 
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The dynamical rule for the ^"■'^ process is that ^"''^ jumps up by 1 at the event times of 
1^i+k provided the inequalities 

er''<ei and C'<ai + 1 (5.32) 
are not violated. In particular ^"''^ attempts a jump at rate c"_,_^. 
Lemma 5.10 (Envelope Property). For each n eN, for all i E Z and t >0, 

zut) = sup{zm - eiAm (5-33) 

feez 

Equation (5.33) holds by construction at time t = 0, and it is proved by induction 
on jumps. For details see Lemma 4.2 in [27]. 

5.2.2 The limit for ^ 

For q,x eM., t > and for the speed function c{x), define 

g'^ix, t) = inf {y : {x, y) G W, r^(a;, y)>t}. (5.34) 

r^(a;,7/) defined by (5.6) represents the macroscopic time it takes a ^-type interface 
process to reach point (x, y). The level curve of given by g'^{-, t) represents the limiting 
interface of a certain ^-process, as stated in the next proposition. 

Proposition 5.11. For all g, a; e M and t > 

lim n-^^";^7J(nt) = ^-'?(x,t) a.s. (5.35) 
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Recall the lattice wedge L defined by (5.1). For (i, j) e £ U let 

L-''=(i,j)=mfO>0:e'W>j} (5.36) 

denote the time when reaches level j. The rules (5.31)-(5.32) give the boundary 
conditions 

L"''=(i,j) = for {i,j)edC (5.37) 
and for G C the recurrence 

L-''{t,j) = max{L-''{t - l,j),L-'\t,j - l),L-''{t + l,j - 1)} + P^f (5.38) 

where P^J^^ is an exponential waiting time, independent of everything else. It represents 
the time ^"''^ waits to jump, after ^^'^ and its neighbors ^^^^ have reached positions 
that permit ^^'^ to jump from j — 1 to j. The dynamical rule that governs the jumps of 
^"l'^ implies that j^lf has rate d^^^. 

Equations (5.4), (5.5), (5.37), and (5.38), together with the strong Markov property, 
imply that 

{L"'^(i,j) : (i,j) e £ua£} ^ {r"'-^(i,j) : {hj) G £ U (5.39) 

Consequently Theorem 5.1 gives the a.s. convergence n~^L"''l-"'«-l([nxJ, \ny\) — >■ r~^(x,y), 
and this passage time limit gives limit (5.35). 

Proof of Theorem 4-2. Given the initial configurations ?7"(0) = {ri^{0) : i G Z} that ap- 
pear in hypothesis (4.12), define initial configurations z"'{0) = {2;"(0) : i e Z} so that 
z'^{0) = so that (5.29) holds at time t^O. Hypothesis (4.12) implies that 

lim n~^zV , = vo(q) a.s. (5.40) 
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for all g e R, with Vo defined by (4.8). 

Construct the height processes and define the exclusion processes r^" by (5.29). 
Define v{x, t) by (4.11). From (5.29)-(5.30) we see that Theorem 4.2 follows from proving 
that for all X e e R+, 

lim n~^z?^^,{nt) ^ v{x,t) a.s. (5-41) 
Rewrite (5.33) with the correct scahng: 

.f,,j(nt) = sup {n-i^f„,j(0) - n-ieH_L.,jM} ■ (5-42) 



lb Z\ 



The proof of (5.41) is now to show that the right-hand side of (5.42) converges to 
the right-hand side of (4.11). 

Prom (5.40), (5.42) and (5.35) we can prove that a.s. 

lim n~^z?^^Ant) = sup {uo(g) - 9~'^{x - q,t)} = v{x,t). (5.43) 



The argument is the same as the one from equations (6.4)-(6.15) in [27] so we will not 
repeat it here. 

Using (5.6) and (5.34) we can rewrite v{x,t) as 

7(x'(.)) 



v{x, t) 



sup (vo(q) -y:3xe n(x - q, y) such that / 7 , j^^^ , ds>t\. (5.44) 
9,2/eM >- 7o c[xi[s) + q) ) 



The final step is to prove v{x,t) = v{x,t). The argument is identical to the one used 
to prove Proposition 4.3 in [28] so we omit it. With this we can consider Theorem 4.2 
proved. □ 
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5.3 Density profiles in two-phase TASEP 

This section proves Corollary 4.4: assuming c{x) = (1 — H(x))ci + H{x)c2, ci > C2 and 
po{x) = p G (0,1), we use variational formula (4.11) to obtain explicit hydrodynamic 
limits. 

Remctrk 5.12. In light of Theorem J^.l, one can (instead of doing the following com- 
putations) guess the candidate solution for the scalar conservation law (4.22) and then 
check that it verifies the entropy conditions (4.24) - (4.26). The following computations 
do not require any knowledge of p.d.e. theory or familiarity with interface problems so 
we present them independently in this section. 



C°{x, t, q) ^ {w e C([0, i], R) : w piecewise linear, w{0) = q, w{t) = x} . (5.45) 

To optimize in (4.11) we use a couple different approaches for different cases. We outhne 
this and omit the details. 

One approach is to separate the choice of the starting point q of the path. By setting 



Let 




(5.46) 



(4.11) becomes 



v{x, t) = sup {vo{q) - I{x, t, q)} . 



(5.47) 



We distinguish four cases according to the signs oi x,q. Set 



R+{x, t) — sup {vo{q) — I{x, t,q)} , if x > 0, 



(5.48) 



q>0 



L_(x, t) — sup {volq) — I{x, t, q)} , ii x < 0. 



(5.49) 



q<0 
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These functions are going to be used in Cases 1 and 2 below {qx > 0) where we can 
compute I{x,t,q) directly. 

However, there are values {x, t, q) for which the g-derivative of the expression in braces 
in (5.47) is a rational function with a quartic polynomial in the numerator. While an 
explicit formula for roots of a quartic exists, the solution is not attractive and it is not 
clear how to pick the right root. Instead we turn the problem into a two-dimensional 
maximization problem. 

li qx < the optimizing path w crosses the origin: w{u) — for some u. It turns 
out convenient to find the optimal q for each crossing time u. For Case 3 {q < 0,x > 0) 
set 

^(u, q)^qp- ciug ( — ) - C2(t - u)g ( ^ ) (5.50) 
\UCiJ \[t-u)c2j 

and 

-L+(x, t)= sup ^{u,q). (5.51) 

g<0,ue[0,t] 

For Case 4 (g > 0, x < 0) the obvious modifications are 

*(u, q)^qp- C2ug { — ] - Ci{t - u)g { ^ ) (5.52) 
\UC2J \{t-u)cij 

and 

R_{x,t) = sup 'if{u,q). 

q>0,ue[0,t] 

Rewrite (5.47) using functions R±, L±: 

v{x,t) = max{R+{x,t), L+{x,t)}l{x > 0} + max{R-{x,t), L_{x,t)}l{x < 0}. (5.53) 

Proof of Corollary 4-4- We compute the functions R±, L±. The density profiles p{x,t) 
are given then by the x-derivative of v{x,t). 
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Case 1: X > 0, q > 0. Since C2 < Ci, the minimizing w of I{x,t,q) is the straight hne 
connecting (0, g) to {t,x). In particular, 



I{x,t,q) = C2tg 



X — q 

tC2 



(5.54) 



Then the resulting R+{x,t) is given by 



-te2^(^) ifp<|, x<tc2{l-2p) 

R+{x,t) = ■lpx_tc2p{l-p), ifp<|, a;>te2(l-2p) 



(5.55) 



px — tc2p{l — p)-, iip> 



2' 



Case 2: X < 0, 5 < 0. The minimizing path w can either be a straight hne from (0, q) 
to (t, x) or a piecewise linear path such that the set {t : w{t) = 0} has positive Lebesgue 
measure. This last statement just says that the path might want to take advantage 
of the low rate at a; = 0. We leave the calculus details to the reader and record the 



4B 



resulting minimum value of I{x, t, q). Set B — -\/ci(ci — C2). 

(i-f)' + (^-^)f-!t(i + |)', 

when - {s/Bi - ^J\x\f < g, -Bt <x<0 
citg (^^^^ otherwise 



I{x,t,q) = < 



(5.56) 
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The corresponding function L_{x,t) is given by 

px — tcip{l — p), < p < e R 

px -tcip{l- p), P*<P<\, x<-tci{p-p*) 

+ + P*<P<h x>-tc,{p-p*) 

px -tcip{l - p), l<p<l-p*, x<-tci{p-p*) 

-(^+l)f + + |<P<1-P*, -tc,{p-p*)<x 

-{t+%)f + m{l + ^y, l-p*<p<l, -Bt<x 

-tcig{^), l-p*<p<l, -cit{2p-l) <x<-Bt 



L_{x,t) = < 



px-citp{l- p), 



1 - p* < p < 1, X < -cit{2p - 1) 



(5.57) 



Case 3: X > 0, q < 0. Abbreviate D — C2 — 4ciC2p(l — p). First compute the g-derivative 



<^>q{u, q) = < 



(5.58) 



p q < —uci. 

If P > 1/2 then $g is positive and the maximum value is when g = so we are reduced 
to Case 1. If p < 1/2 the maximizing q — uci (2p — 1) . Then 

F{u) = <^{u, 2uci{p - \)) = -wcip(l - p) - C2{t - u)g -u)o^ ' 



with li-derivative 
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Again we need to split two cases. If p < p* (equivalently D > 0) and x < WD, the 
maximizing u — t — x/\/D, otherwise u — 0. li p* < p < ^ the derivative is negative so 
the maximizing u is still u — 0. Together, 



L+{x,t) = < 



-tcip(l-p)+x(|-|f), p<p*,x<tVD 



p < p*,x > tVD 



P* < P<1- 



(5.59) 



Case 4: X < 0,q > 0. We treat this case in exactly the same way as Case 3, so we omit 
the details. Here we need the quantity Di — (ci)^ — 4ciC2p(l — p) and we compute 



X 



(5.60) 



R-{x, t) = { _tc^p(i _ p) + ^ 1 + ^ j , 1 < _tyfDl < 
Now compute v{x,t) from (5.53). We leave the remaining details to the reader. □ 



5.4 Entropy solutions of the discontinuous conser- 
vation law 

For this section, c{x) — {1 — H(x))ci+H(x)c2, h{p) — p{l—p) and set F{x,p) — c{x)h{p) 
for the flux function of the scalar conservation law (4.19) and F{x, p) = c{x)f{p) for the 
flux function of the particle system, where / is given by (4.9). (The difference between 
F and F is that the latter is — oo outside < p < 1.) 
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In [1] the authors prove that there exists a solution to the corresponding Hamilton- 
Jacobi equation 

Vt + Cih{V^) =0, if X < 0, t > 
' Vt + C2h{V^)=Q, ifx>0,i>0 (5.61) 
^ Vix,0)^Vo{x) 

such that Vx solves the scalar conservation law (4.19) with flux function F{x,p) and 14 
satisfies the entropy assumptions {Ei), {Ef,). V{x,t) is given by 

V{x,t)^sup{vo{w{0))+ [ {c{w{s))hy{w'{s))ds \ , (5.62) 
w{-) [Jo ) 

where the supremum is taken over piecewise linear paths w e C([0,t],M) that satisfy 
w{t) — X. 

To apply the results of [1] to the profile p(x, t) coming from our hydrodynamic limit, 
we only need to show that the variational descriptions match, in other words that we 
can replace F with F and the solution is still the same. 

Proof of Theorem 4-7. Convex duality gives {c{x)f)* (y) — c{x) f* {y / c{x)) and so we 
can rewrite (4.11) as 

v{x, t) = sup lvo{w{0)) + / {c{w{s))f )*{w'{s)) ds\ . (5.63) 
w{-) { Jo ) 

Observe that for all y e R 

{c{x)fr{y) > {c{x)hny), (5.64) 

with equality if and only if y G [— Ci, C2] Since the supremum in (5.62) and (5.63) is taken 
over the same set of paths, (5.64) implies that 

V{x,t) <v{x,t). (5.65) 
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The proof of the theorem is now reduced to proving that the supremum in (5.63) is 
achieved when w'{s)c{w{s))~^ G [—1, 1], giving V{x,t) = v{x,t). 
To this end we rewrite v{x,t) once more, this time as 

v{x, t) — max{i?+(a;, t), L+{x, t)}l{x > 0} + max{R-{x, t), L_(x, t)}l{x < 0} 

where the functions R±, L± (as in the proof of Corollary 4.4) are defined by 

R+{x, t) — sup {vo{q) — I{x, t, q)} , if x > 0, (5.66) 

q>0 

L_{x, t) = sup {vo{q) — I{x, t, q)} , if x < 0, (5.67) 
where I{x,t,q) is as in (5.46), and 

L+{x,t) ^ sup \vo{q) - ciug(—) - C2{t - u)g( ^ )\ if x > 0, (5.68) 

q<0,u&[0,t] I ^UCiy \[t-U)C2y } 

and 

R^{x,t)^ sup \vQ{q)-C2ug( — ]-Ci{t-u)g( ^ ] \, x<0. (5.69) 
g>0,uelO,t] [ \uc2J \{t-u)cij } 

It suffices to show that the suprema that define R±, L± are achieved when 

w'(s)c(w(s))-^ e [-1,1]. (5.70) 

We show this for L+. The remaining cases are similar. In (5.68), as before, u is the 
time for which w{u) = 0. Let $(«, q) denote the expression in braces in (5.68) with 
g-derivative 

Poi<l) - i - 2^' <q<0 
%{u,q)^l (5.71) 

Po{q), q < -uci. 
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Observe that if ^q{u, q) — for some q* — q*{u) then also q* maximizes Otherwise 
the maximum is achieved at and we are reduced to a different case. Assume that q* 
exists. Then by (5.71) 

^ = (l-2po(9*))ci<Ci. (5.72) 
u 

Therefore, the slope of the first segment of the maximizing path w satisfies (5.70). 

The slope of the second segment is x{t — u)~^. Assume that the piecewise hnear path 
w defined by u and q* is the one that achieves the supremum. Also assume u > t — xc^^. 
Consider the path w with ^(0) = q* and w{t — xc^^) — 0. Since g is decreasing, we only 
increase the value of $. Hence the supremum that gives L_|_ cannot be achieved on w 
and this gives the desired contradiction. □ 
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Appendix A 



Basic Facts 



In this section we report all basics facts from analysis, special functions and probability 
theory used throughout the thesis. 



We only use it for positive real values of s. 

The logarithm logr(s) is convex and infinitely differentiable on (0, oo). The deriva- 
tives are called polygamma functions ^„(s) = {d"''^^ / ds^''^^) logr(s), defined for n G Z+. 
For n > 1, ^„ is nonzero and has sign (— l)"'"^ throughout (0, oo). In particular, ^o{s) 
is strictly increasing and has a vertical asymptote at s = 0. It can be given by 



One way to compute the hmit (3.19) is by multiple uses of L' Hospital's rule and 
then an asymptotic analysis for ^i(s) for s ^ 0. For the asymptotic analysis, we need 



A.l Special functions and distributions 



The gamma function is 




oo 




(A.l) 




Ill 

The Gamma(^, 1) distribution has density 

r(e)-^x^-^e-^, e>0. (A.2) 

on R+, mean 9 and variance 9. 

Throughout the dissertation we make use of the digamma and trigamma functions 
^0, ^1 since for X ~ Gamma(^, 1) we have 

*o(^) =E(logX) and = Var (log X). (A.3) 

A.2 Convex Analysis 

For given functions f{x),g{x) we denote the convex dual 

fir) = sup{rx - fix)}, (A.4) 

X 

and the infimal convolution 

{fDg){x) = inf{/(y) + g{x - y)}. (A.5) 
y 

For lower semi-continuous convex / and g we have 

{fngr = r+g*, (A.6) 

and double convex duality 

/** = /■ (A.7) 
Also, convexity of / implies that on the set {|/| < oo}, / is a.e. differentiable with 

f{x) = axgmax{xu - r{u)}. (A.8) 
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A. 3 Large deviations 

Here are some basic theorems from the theory of large deviations that we are using 
throughout. The hmiting log-moonet generating function is given by 

M(m) = hm n-4ogE(e"^-i^^) (A.9) 

Theorem A.l (Cramer's Theorem). Let {Xn}n be a sequence of i.i.d. real-valued ran- 
dom variables. Let /in be the law of the sample mean Sn/n. Then, the large deviation 
principle LDPliinjn, I) is satisfied with I defined by 

I (a) = sup{aM - M{u)}, (A. 10) 

where M{u) is the limiting log-moment generating function given by (A.9). 

Theorem A. 2 (One sided Cramer's Theorem). Let {Xn}n be a sequence of i.i.d. real- 
valued random variables. Define the one sided rate functions by 

J{a) = hm n~MogP{,S„ > na), (A.ll) 



/(a) = hm n~^\og¥{Sn < na}. (A. 12) 

n->oo 

The the two functions are given respectively by 

J{a) = sup{aw - M(w)}, (A.13) 

u>0 

I{a) = sup{au- M{u)}. (A.14) 

where M{u) is given by (A.9) 
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